Unit 8 


Differential equations 


Introduction 


Introduction 


This unit introduces differential equations. A differential equation is an 
equation that includes derivatives — here are two examples: 
2 
cu = 2y and Z + ort + 9a7y = sin(2z). 
In both these equations y is being differentiated with respect to x, so y is a 
function of x. The process of finding y as a function of x is called ‘solving 
the differential equation’. 


Differential equations arise in many situations in mathematics and the 
mathematical sciences, such as physics, chemistry, engineering, biology and 
even economics and finance. This is because the rates of change of the 
quantities involved (which are derivatives) often behave in much simpler 
ways than the quantities themselves. Some important applications of 
differential equations are mentioned at the end of this introduction. 


This unit concentrates on three types of differential equations, the simplest 
one of which is a special case of each of the other two. After studying the 
unit, you should be able to recognise differential equations of each of the 
three types, and apply an appropriate method to solve each type. 


Section 1 deals with differential equations of the form 


Such equations can be solved by direct integration. 


Section 2 shows how differential equations that can be written in the form 


dy 
ag 7 fa) 
can be solved by using the method of separation of variables. 


Section 3 investigates a number of applications of differential equations to 
real-world problems. It differs from the other sections in this unit, in that 
here you will see how differential equations can arise in applications. The 
section goes on to solve these differential equations and explore some 
properties of their solutions. 


Section 4 deals with differential equations of the form 

d 

T + gley = h(a). 

These equations are called linear differential equations, and are solved 


using a ‘trick’ that turns them into a form that can be solved by direct 
integration. 


Section 5 summarises and revises all the methods for solving differential 
equations that are covered in the unit. Finally, in Section 6 you’ll see how 
differential equations can be solved on a computer. 
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Differential equations in science 


In physics, differential equations often arise from fundamental 
physical laws. Here are four important examples. 


e Newton’s second law of motion is essentially a differential 
equation that describes how objects move under the influence of 
forces. 


e Mazwell’s equations are differential equations that describe the 
behaviour of electricity and magnetism. 


e Finstein’s equations of general relativity are differential equations 
that describe the motion of extremely massive objects interacting 
Isaac Newton (1642-1727) via gravity, such as stars, galaxies and black holes. 


e Schrödinger’s equation is a differential equation that describes 
the motion of particles on the atomic and subatomic scale. 


Other fundamental physical laws are also described by differential 
equations, and together these differential equations accurately 
describe the physics and chemistry of the vast majority of the 
Universe as we know it. 


Differential equations also arise in subjects such as biology, 
engineering and economics, and also in physics and chemistry, to 
describe the behaviour of very large numbers of individual objects in 
ageregate, when the complexity of describing every individual in the 
group is too great. 


James Clerk Maxwell For example, the motion of water could be described in principle by 
(1831-1879) describing the motion of each individual water molecule. However, 
since there are approximately 3 x 10? molecules in a single gram of 
water, the resulting system of differential equations would be 
extremely complicated. Instead, to a very good approximation, the 
motion can be described by a small number of differential equations, 
together called the Navier—Stokes equations. 


Other phenomena that can be described in this way include: the 
weather, the spread of disease, population changes, electrical circuits 
and the economies of countries. 


1 Differential equations and direct 
integration 


This section introduces differential equations and the idea of solving a 
differential equation. It then illustrates these concepts by looking at the 


j p simplest type of differential equation, namely those that can be solved by 
Erwin Schrödinger direct integration. 
(1887-1961) 
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1 Differential equations and direct integration 


1.1 Differential equations 


In general, a differential equation is an equation involving an unknown 
function and its derivatives. In the introduction to this unit you saw the 
following two examples: 


dy 

== 9 

dx Ys (1) 
and 

q? d 

T y% 6x + 9z?y = sin(2z). (2) 


In each of these differential equations the variable y is a function of 2; 
we sometimes denote this fact by writing y = y(x). Since y is a function 
of x in these equations, y is the dependent variable and x is the 
independent variable. Note that the dependent variable y is 
differentiated with respect to the independent variable x. 


Letters other than x and y are sometimes used for the independent and 
dependent variables in a differential equation, especially t for the 
independent variable to represent time. For example, 


d d 

“4 _ #2 is the same differential equation as peasy : 

dt da 
except that the independent variable has changed from x to t and the 
dependent variable from y to q. 


2 


The nth derivative of a function is said to be the derivative of order n, 
and the order of a differential equation is the order of the highest 
derivative that appears in the equation. For example, differential 
equation (1) is of first order since the highest derivative that appears is 
dy/dx, and differential equation (2) is of second order since the highest 
derivative that appears is d?y/dz?. In this unit you will meet only 
first-order differential equations. More specifically, you will meet only 
first-order differential equations that can be expressed in the form 

% = fley), 
where f(x,y) denotes an expression that depends on the variables x and y. 
For example, the differential equations 


dy _ dy _ 2 
dn otY and ie +y°, 


are both of this form. 


A solution of a differential equation is a function y = y(x) that satisfies 
the differential equation when substituted into that equation. For example, 


the function y = e?” is a solution of the first-order differential equation 
dy 
— = 2y, 3 
da Y (3) 
because with y = e?” the left-hand side of the differential equation equals 
dy d / 2¢ 2 
ea een 


which is the same as the right-hand side. 
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Usually a differential equation has more than one solution. For example, 
two other solutions of differential equation (3) are the functions y = 3e?” 
and y = —5e?*, as you can easily check. In fact, any function of the form 
y = Ae?*, where A is a constant, is a solution of equation (3). 


To solve a differential equation you usually try to find the complete 
family of solutions, expressed in the form of a general solution involving 
one or more arbitrary constants. For example, as you will see later in the 
unit, the general solution of equation (3) is y = Ae?*, where A is an 
arbitrary constant. 


If you choose a particular value for each arbitrary constant in the general 
solution of a differential equation, then you obtain a particular solution 
of the differential equation. For example, y = e°”, y = 3e?” and y = — 5e?” 
are all particular solutions of equation (3). 


1.2 Direct integration 


The simplest type of differential equation is one of the form 


Here the right-hand side is an expression in the independent variable x and 
contains no terms involving the dependent variable y. We call differential 
equations of this type directly integrable. 


For example, the differential equations 


Le dy Sag? Uys To 
dx” dx ~? dx i dr z 
are all of this type, but 
w = ry and w =x“-y 
are not. 


Activity 1 /dentifying directly integrable differential equations 


Which of the following differential equations are directly integrable? 
da 


dy dx 
(a) qn T Ecos (b) qe T7082 (c) qg 70st 


You can solve a directly integrable differential equation 


= = f(x) (4) 


simply by integrating both sides with respect to the independent variable, 
provided that you can integrate f(x). 


1 Differential equations and direct integration 


Since integration reverses differentiation, this gives 


y= | fæða, 


so the solution of equation (4) is just the indefinite integral of the 
function f. This method is called direct integration. 


Let’s apply direct integration to a simple example. 


When you have solved a differential equation, it’s often a good idea to 
check that your solution does satisfy the original differential equation. 


With y = $x + c in Example 1, the left-hand side of the differential 
equation is 

dy d 13 2 

o +c) =2", 
which is the same as the right-hand side, so the differential equation is 
satisfied. 


Here’s another example. 
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Example 2 Solving another differential equation by direct 
integration 
(a) Solve the differential equation 


dy 3 


(b) Check that your solution satisfies the differential equation. 


Solution 


(a) @ The differential equation has the form dy/dx = f(x), so it can 
be solved by direct integration. ® 


Integrating both sides gives 
y= [@ -542 = let _5rte, 


where c is an arbitrary constant. 


(b) With y = 4x1 — 5x + c the left-hand side of the differential 
equation is 
dy drm 3 
which is the same as the right-hand side, so the differential 
equation is satisfied. 


Here’s an example of direct integration for you to try. 


Activity 2 Solving a differential equation by direct integration 


(a) Use direct integration to solve the differential equation 


dy 
— = 2 3. 
Tr x + 


(b) Check that your solution satisfies the differential equation. 


General solutions obtained by using direct integration always contain an 
arbitrary constant. For example, in Example 1 you saw that the 
differential equation 


dy 2 

an (5) 
has the general solution 

y = že’ + C, (6) 


where c is an arbitrary constant. If you set the arbitrary constant in a 
general solution to be a particular value, then you obtain a particular 


1 Differential equations and direct integration 


solution. For example, if you set c = 5 in equation (6), then you obtain the 
function y = ir’ + 5, which is a particular solution of differential 
equation (5). 


Here is an example where you are asked to identify particular solutions. 


Example 3 Finding a general solution and particular solutions 
(a) Find the general solution of the differential equation 
= = ¢7". 
(b) Which of the following functions are particular solutions of this 
differential equation? 
(i) y= ge +7 (ii) y= se — (ili). y = Ge? +30 
(iv) y = 5(e?* + 3) 
Solution 


(a) ®. The differential equation has the form dy/dz = f(x), so it can 
be solved by direct integration. © 


Integrating both sides of dy/dx = e?” gives the general solution 
y= {2 Gre ek, a= ze" +c, 


where c is an arbitrary constant. 


@. Check: differentiating y = Je?” + gives dy/dx = e?”, as 
required. ® 


(b) The functions in (i), (ii) and (iv) are particular solutions obtained 
by setting the constant c equal to 7, 0 and 3, respectively. The 
function in (iii) is not a particular solution as it is not of the form 


y = ġe” +c, where c is a constant. 


Here is a similar example, in which t is the independent variable. 


Activity 3 Finding a general solution and particular solutions 


d 
(a) Find the general solution of the differential equation = = sind. 
(b) Which of the following functions are particular solutions of this 
differential equation? 


(i) c=2cost+1 (ii)z=cost—3 (iii) z= 1-— cost 
(iv) z = — cost + 2t 
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When using differential equations, you often need to find a particular 
solution that has a given additional property. To illustrate the ideas, 
consider once again the differential equation from Example 1, 

d 

dy _ 42 

dx 
which has the general solution 

y= r+, 
where c is an arbitrary constant. Different values for c give different 
particular solutions, and graphs of the particular solutions corresponding 


toc= $, z, 1 and - are shown in Figure 1. 
e= g 
G= 
C= 23 
e= 1/3 


SY 


| 
| 
-1 0 1 


Figure 1 The graphs of y = ir’ +c, for c = $, z, 1 and $ 


Now, you might want to choose the particular solution whose graph passes 
through the point (1,1). As shown in Figure 1, the particular solution with 
this property is 

y= ix’ + 2, 
Notice that if you choose any point in the plane, then only one of the 
graphs of the form y = ir’ + c passes through this point. So choosing a 


point through which the graph of the solution must pass is equivalent to 
picking a particular solution. 


Requiring that the solution to a differential equation passes through a 
given point in the plane is a very common way of specifying a particular 
solution. If a particular solution is specified in this way, by requiring that 
its graph passes through a given point, then this requirement is called an 
initial condition. In the example above, the particular solution 

y= ig? + 3 satisfies the initial condition y = 1 when x = 1. 

Here is an example of a problem that uses an initial condition to specify a 
particular solution. 
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Example 4 Applying an initial condition 

In Example 2 you saw that the differential equation 
“ =f —5 

has the general solution 
y= Izt = Sue e 

where c is an arbitrary constant. 


Find the particular solution of this differential equation that satisfies 
the initial condition y = —1 when x = 2. 


Solution 


@. We need to find the specific value of c such that y = —1 when 
c=2. ®@ 


We substitute the values x = 2 and y = —1 (simultaneously) into the 
equation y = 424 — 5a +c. This gives 


-1l=4x2*-5x2+c=-6+¢, 
so c= 5. The required particular solution is therefore 


y= lr -5r +5. 


An initial condition can be written using the y(x) notation for the 
unknown function. For example, you can write the initial condition 


y = —1 when z = 2 
more simply as 


y(2) =-1. 


Here is a similar problem for you to try. 


Activity 4 Applying an initial condition 


The differential equation 
dy 5 


dr 2./x 


has general solution 


y=5/a+¢, 


where c is an arbitrary constant. Find the particular solution of this 
differential equation that satisfies the initial condition y = 7 when x = 1; 
that is, y(1) = 7. 
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Here is a summary of the method of solution by direct integration. 


Direct integration 


The general solution of the differential equation 


dy _ 
ap = (x) 


is the indefinite integral 


y= | f@)ax, 
which can be written as 
y=F(z)+¢, 


where F(x) is any antiderivative of f(x) and c is an arbitrary 
constant. 


An initial condition of the form 
y=bwhenz=a, thatis, y(a)=68, 


enables you to find a value for the arbitrary constant c. The 
corresponding particular solution satisfies both the differential 
equation and the initial condition. 


This method can always be applied to solve a differential equation of the 
form dy/dx = f(x), as long as you can find the indefinite integral involved. 


A problem in which you have to find the particular solution of a 
differential equation that satisfies a given initial condition is called an 
initial value problem. Generally there are two stages in solving an 
initial value problem. 


Solving an initial value problem 
To solve an initial value problem: 
1. Find the general solution of the differential equation. 


2. Use the general solution to find the particular solution that 
satisfies the initial condition. 


Here is an example of solving an initial value problem. In this example the 
differential equation has a restriction t > 0. This restriction is given since 
the expression 1/t? is not defined when t = 0, so we look for a solution 
restricted to an interval that contains the value of t in the initial condition. 
This is because when we are solving a differential equation we are usually 
interested only in solutions that are continuous functions. 
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Example 5 Solving an initial value problem 


Solve the initial value problem 


d 1 
=F ({>0), where g=1 when t= 2. 
Solution 


@. First find the general solution. The differential equation has the 
form dq/dt = f(t), so it can be solved by direct integration. ® 


Integrating both sides of the differential equation with respect to t 
gives the general solution 


1 
q= ad= frra=- +e (t > 0), 


where c is an arbitrary constant. 
@®, Substitute the initial condition to find c. ® 


Since q = 1 when t = 2, we obtain 
1 : F 
=--+¢ 
9) ? 


which gives c = 3. Hence the particular solution that satisfies the 
initial condition is 

— l 3 0 

¢G = 7 T 3 (t > ye 


Here are some similar examples for you to try. 


Activity 5 Solving initial value problems 


Solve each of the following initial value problems. 
dy 


(a) as yx, where y=5 when z = 4 
x 
du 1 
(b) Ti cos(2x), where u = —2 when x = 57 


Sometimes a whole family of similar differential equations can be solved at 
the same time by including a letter that denotes a constant in the 
equation; such a letter is often called a parameter. Here’s an example of 
a differential equation with a parameter for you to try to solve. 


© 


165 


Unit 8 Differential equations 


A model Saturn V rocket on 
launch 
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Activity 6 Solving an initial value problem with a parameter 
Solve the initial value problem 


d k 
= =o (t >0), where k is a positive constant, and x = 3 when t = 1. 


The final activity in this section involves a possible application of 
differential equations. In this application the derivative dm/dt represents 
the rate of change of the mass, m, of a rocket with respect to time, t. The 
mass of the rocket is decreasing with time as the fuel is used up, and hence 
its rate of change is always negative. 


Activity 7 An application: finding the mass of a burning rocket 


A toy rocket consists of a casing and fuel. Initially the total mass of the 
rocket is 200 grams. The fuel is ignited at time t = 0 and then burns in 
such a way that the rocket’s total mass m (in grams) at time t (in seconds) 
satisfies the differential equation 


dm 


a > 0). 
= =t (t20) 


(a) Find the general solution of the given differential equation. 


(b) Find the particular solution of the differential equation that satisfies 
the initial condition m = 200 when t = 0. 


(c) When the fuel is exhausted, the residual mass of the rocket casing is 
100 grams. Find the time at which the fuel is exhausted; that is, find 
the time t (in seconds) when m = 100 for the particular solution 
obtained in part (b). 


2 Separable differential equations 


In this section we consider a second type of first-order differential equation, 
called a separable differential equation, and we show how to solve separable 
differential equations to find the general and particular solutions. 


2.1 Separation of variables 


Direct integration applies only to the very simplest type of differential 
equation, those of the form dy/dx = f(x). In this subsection, we consider 
how to solve first-order differential equations of the form 


Y = F(a)gly), (7 


2 Separable differential equations 


where the right-hand side is the product of a function of x and a function 
of y. We can’t solve such equations by direct integration, because the 
right-hand side depends on both x and the unknown function y = y(x). 


Differential equations of the form in equation (7) are called separable. 
For example, 


dy 2 y . dy zx dy siny 

— =f — = sin 7 cosy, — = >, a 

dx 2 dx Y dz y? dz  cosz 
are all separable, but 

dy dy 2 dy 

dx TY dc” Tai dz n(ry), 


are not separable because their right-hand sides cannot be written as a 
function of x multiplied by a function of y. 


Activity 8 /dentifying separable equations 


Identify which of the following differential equations are separable; that is, 
which of them are of the form 
d(dependent variable) 


deaa] = f(independent variable) g(dependent variable). 


For each equation that is separable, identify the corresponding functions f 


and g. 

(a) Gp T 7CSy (b) g T cost (c) qe PoP 
dv dv | dy 

Gg rte pata Op te 


To illustrate the method for solving a separable differential equation, let’s 
start with a fairly simple example. Consider the differential equation 

dy 2? 

=Z (y>0). (8) 

T y 

This is of the form in equation (7), with f(a) = x? and g(y) = 1/y, so it is 
a separable differential equation. The condition y > 0 is included because 
the right-hand side of the differential equation is not defined when y = 0. 


To solve differential equation (8), first rearrange it so that the dependent 
variable and its derivative appear only on the left-hand side and the 
independent variable appears only on the right-hand side. Multiplying 
both sides of the equation by y gives 


dy > 


Now integrate both sides with respect to x: 


fo Zas fean (9) 
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Next recall the rule for integration by substitution from MST124, revised 
in Section 1 of Unit 7, which tells you that any integral of the form 


frost dz (10) 


is equal to the simpler integral 


f Hia 


(It’s easy to remember this rule because it looks as if we simply cancel a 
‘dx’ in the denominator with a ‘dx’ in the numerator, though of course 
that’s not what’s actually happening.) 


By applying this rule to the left-hand side of equation (9), which is of the 
same form as expression (10) with y instead of u, we can simplify 
equation (9) to 


[veu= [oar (11) 


Now you can see that the variables x and y have been ‘separated’, so you 
can integrate both sides of this equation to give an equation relating y to x. 


Doing the integrations on both sides gives 
ay” = ig’ FG 

where c is an arbitrary constant. 

(You could have introduced an arbitrary constant on both sides to obtain 
sy +a = 30° Ho, 

but the constants cı and c2 could then be combined to give iy? = ir? +c, 

where c = c2 — c1.) 

Hence 


y = 44/323 +C, (12) 


where C = 2c is another arbitrary constant. 


Now remember the condition y > 0 that applied to the original differential 
equation. This tells you that the solution y = y(x) must take only positive 
values. Hence the general solution of the differential equation, with this 
condition, is 


y=/30+C. (13) 


At this point you could check this solution by differentiating it and 
substituting it back into the original differential equation, as in 
Subsection 1.2. This would be slightly laborious but perfectly doable. 
However, later in this section you will learn a technique for checking your 
answer that is usually more straightforward. 


2 Separable differential equations 


Notice that the general solution found above contains an arbitrary 
constant, as expected. However, unlike the general solution of a directly 
integrable differential equation, this general solution is not of the form 

y = F(x) +c, since the arbitrary constant in it cannot be taken outside the 
square root. 


You can find particular solutions of separable differential equations in 
exactly the same way as for directly integrable differential equations. 
Suppose, for example, that you want to find the particular solution of 
equation (8) that satisfies the initial condition 


y = 2 when z = 3. 


To do this you simply substitute x = 3 and y = 2 into the general solution, 
which gives 


2= 1/2 x33 +C. 


Squaring both sides gives 
2? = 2x3 +O, 


so C = —14. Hence the required particular solution is 


y = 4/323 — 14. 


Notice that this particular solution is a function that is not defined for all 
values of x. We shall return to this point later in the section. 


Before you see some further examples of separating variables, here is a 
shortcut that makes the process easier. In the discussion above, we used 
integration by substitution to show that if 


dy x? 


dr y’ 
then 


[ucy= far 


A convenient way to remember this step is to imagine that the derivative 
dy/dz is a fraction, at least temporarily. Then you separate the variables 
by taking everything involving y to the left-hand side of the differential 
equation and everything involving x to the right-hand side, and finally 
putting integral signs on both sides. 


This shortcut works for any differential equation of the form in 
equation (7) and is summarised in the following box. 
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A shortcut for separating variables 


To get from the differential equation 
dy 
Z = se) lv) 


to 
[gr] o 


imagine that dy/dz is a fraction, and take everything involving y to 
the left and everything involving x to the right, and finally put 
integral signs on both sides. 


The Y's to the Left... 
The Y's have it! 


You should remember however, that the derivative dy/dz is not really a 
fraction, but the clever notation dy/dx makes this manipulation appear 
natural. 


Once you have separated variables in this way, various possibilities can 
occur when you try to proceed. 


e You may be able to find both integrals and rearrange the resulting 
equation to express y explicitly in terms of x. The resulting equation 
is a solution in explicit form. 


e You may be able to find both of the integrals, but the resulting 
equation relating y and x may be so complicated that you cannot 
rearrange it to express y explicitly in terms of x. In this case, the 
equation relating y and is said to be a solution in implicit form. 


e One or both of the two integrals may be difficult or impossible to find. 
In this case the method has failed. 


Here is an example that demonstrates the use of the shortcut for 
separating variables and also illustrates the second possibility above. 
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Figure 2 shows the curves in the x, y-plane given by the equation 

y +lny =—cosx+c, for c = —1, 0 and 1. Each of these curves seems to 
be the graph of a function since each value of x appears to correspond to 
just one value of y. This is in fact the case, but we shall not investigate 
this further in this unit. 


=o = a ae a a” 


Figure 2 The graphs of y+Iny = — cos x + c, for c= —1, 0 and 1 


In seeking to solve a differential equation, you should aim to find a solution 
in explicit form wherever possible. However, if this is impossible, as in 
Example 6, then you must settle for a solution in implicit form. 


If you have found the general solution only in implicit form, then you can 
still find a particular solution (in implicit form) from a given initial 
condition. Here’s an example. 
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Example 7 Finding a particular solution in implicit form 
In Example 6 you saw that the differential equation 
dy — ysinz 


da Wap ll 


has general solution 


y +lny =—cosx+c, 
in implicit form, where c is an arbitrary constant. 
Find the particular solution that satisfies the initial condition 


vis event = 0k 


Solution 


Substituting the initial condition y = 1 when x = 0 into the general 
solution y + ln y = — cos x + c, we obtain 


mle e 


so c = 2 + ln 1 = 2. Hence, the particular solution in implicit form 
that satisfies the given initial condition is 


y +lny = —cosz +2. 


The method for solving differential equations that you have met in this 
section is known as separation of variables and is summarised below. 


Separation of variables 


This method applies to differential equations of the form 
dy 
eS . 16 
qa 7 f(@) oy) (16) 
1. Separate the variables and put integral signs on both sides to 
obtain 


[iv [fox 


2. Carry out the two integrations, introducing one arbitrary 
constant, to obtain the general solution in implicit form. 


3. If possible, manipulate the resulting equation to make y the 
subject, thus expressing the general solution in explicit form. 


In some cases it is not possible to carry out step 2 or step 3. 


2 Separable differential equations 


Once a general solution has been obtained, in implicit form or explicit 


form, it can be checked by substituting back into the differential equation. 


You will see how to do this for solutions in implicit form in the next 
subsection. 


Note that if g(y) = 1 in equation (16), then the separable differential 
equation becomes 


oH = f(e). 


In this case the method of separation of variables reduces to simple direct 
integration, the method covered in Subsection 1.2. Hence direct 
integration is just a special case of separation of variables. 


Here are some problems involving separation of variables for you to try. 


Activity 9 Finding a solution in explicit form 


(a) Find the general solution in explicit form of the differential equation 


dz y 
(b) Find the particular solution of this differential equation that satisfies 
the initial condition y = 1 when x = 0. 


Activity 10 Finding a solution in implicit form 


(a) Find the general solution in implicit form of the differential equation 
dy 1+e* 
= = ——__ (-1<y<7). 
dx 1+ cosy nt) 


(b) Find the particular solution of this differential equation that satisfies 
the initial condition y(0) = 0. 


2.2 Verifying solutions in implicit form 


In Subsection 1.2 you saw how to verify that an explicit solution obtained 
by direct integration does satisfy the original differential equation. You 
will now see how to check that a solution you have found in implicit form 
satisfies the differential equation from which it was obtained. 


To do this you'll need to use the chain rule for differentiation, which you 
revised in Unit 1. 
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Chain rule 
If y is a function of u and u is a function of x, then 


ay _ dy du 
dx du daz’ 


The chain rule is needed for checking solutions to differential equations 
when you have found a solution in implicit form involving the independent 
variable x and the dependent variable y, and you want to check that this 
solution satisfies the original differential equation, which involves dy/dx. 
In that situation, you may need to differentiate a function of y such as 


G(y) =y’, with respect to a, 
so 
G is a function of y and y is a function of x. 


In this case the chain rule gives 
dG _ dG dy _ yd 
dz dy dz da 

or, equivalently, 


dpa, dpa dy _ dy 
ae ag ae ae 


Here is an example. 
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By the chain rule, the left-hand side is 
d dys d 
a ee) ar ae”) 
dy dy 
ey a RA 
=e) dz a dz 


The right-hand side is 


Rearranging gives 


dy 
Ta Dye 
that is, 
dy č 2z 
de 3y2 +1’ 


which is the required differential equation. 


Here are two similar problems for you to try; recall that in Activity 10, you 
were asked to solve the differential equation in Activity 11. 


Activity 11 Verifying a solution given in implicit form 


Show that 
ytsiny=x+e"—-1, 
is a solution in implicit form of the differential equation 


dy 1 -pe* 


= —— =—7 YY < 7): 
dz 1+cosy ( d ) 


Activity 12 Verifying another solution given in implicit form 


Show that 
y + lny = (1 — 22)e** 

is a solution in implicit form of the differential equation 
dy — y(1—6a)e** 


de" y+1 «(> 
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A note of caution 


Some care is required when applying separation of variables. First, you can 
divide both sides of the equation 


Y L Fa) aly) 


by g(y) only if g(y) 4 0, since division by zero is undefined. Therefore, 
when separating variables, it is necessary to exclude from consideration 
any value of y for which g(y) = 0. This can be done by imposing a suitable 
condition on y from the outset. For example, in Activity 9 we applied the 
condition y > 0 to the differential equation dy/dx = —a/y, and obtained 
the positive general solution y = vb — z?. As you saw earlier, conditions 
like this often arise from the context of the problem being considered. 


In addition to restrictions on the dependent variable, sometimes solutions 
may be valid only for a certain range of values of the independent variable. 
For example, in Activity 9(b) you found the particular solution 

y = V1 — x?. For this function to be defined we require 1 — x? > 0, which 
is equivalent to —1 < x < 1. Since x = 1 and x = —1 both give y = 0, 
which is excluded by the restriction y > 0 in the question, the domain of 
this particular solution is —1 < x < 1. This example shows that 
consideration of the appropriate domain for a solution can sometimes be a 
bit complicated. 


Activity 13 Finding a solution in explicit form of an initial value 
problem 


Solve each of the following initial value problems, giving your answers in 
explicit form. 


er 


(a) 7 —y? (y>0), where y= 5 when z = 0. 
z 
u? du 
(b) Aao (u >0, x>0), where u = 1 when z= 1. 


Activity 14 Finding a solution in implicit form of an initial value 
problem 


Solve the following initial value problem, giving your answer in implicit 
form. 
dy _ y? cosg 
dx yt+l 


(y>0), where y = 2 when z = z 


2 Separable differential equations 


In the final activity of this section you will use separation of variables to 
solve a practical problem concerning water flow, in which the differential 
equation is derived from a physical law called Torricelli’s theorem. In this 
activity the differential equation involves a parameter k and the initial 
condition contains a parameter yo. Remember that k and yo are constants. 
The notation yo is used to denote the value of the variable y at the initial 
time t = 0. 


This activity concerns a mathematical model for the flow of water from a 
tank, of constant horizontal cross-section, through a small hole at its base. 
This is shown in the figure below. 


According to the model, the height y of the water surface above the hole is 
related to the time t that has elapsed since the hole was uncovered. This 
relationship is described by the differential equation 


dy 

= = ky"? (y>0), 

ET y= (y>0) 

where k is a positive constant. 

(The condition y > 0 is natural in the context of the problem and the 


value of k depends on the area of the hole and the horizontal 
cross-sectional area of the tank.) 


(a) Find the general solution in implicit form of this differential equation. 


(b) Find the particular solution for which y = yo at t = 0. (Here yo is the 
height of the water column above the hole at the moment when the 
hole is uncovered.) 


(c) Show that, according to this model, the tank becomes half-empty after 
(2 — v2) Vm 
—; 


time 
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Figure 3 The world 
population from 1800 to 2010 
(for the years up to 1950, 
based on US Census Bureau 
historical estimates) 
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3 Applications of differential equations 


In this section we will investigate some applications of differential 
equations to real-world problems. You will see how differential equations 
arise in a variety of contexts, and go on to solve these differential equations 
and explore some properties of their solutions. 


We will concentrate on real-world problems that can be modelled by 
first-order differential equations of the form 


dy 
SER 
T y +C, 


where K and C are constants. A differential equation of this type can be 
solved using separation of variables, and the resulting behaviour is that of 
exponential growth or decay, as you will see. 


The variety of applications of differential equations that you will see 
indicates how ubiquitous they are throughout the mathematical sciences. 


3.1 Continuous model for population change 


In this subsection you will see how to use a differential equation to model 
the variation in the size of a population over time. The population could 
be the people in a country or in some other geographical region, or it could 
be the bacteria in a Petri dish, or any other population of individual 
organisms. 


Since a population is made up of individuals, the size of a population can 
take only integer values, so the graph of the population size against time 
will not be continuous. However, if the population size is large, then its 
behaviour over time can be approximated by a function with a continuous 
graph — in other words, a continuous function. For example, Figure 3 
shows the graph of a continuous function that approximates how the size 
of the population of the world has varied with time. 


If you approximate a population size over time by a continuous function, 
then you can consider the rate of change of the population size at any 
particular point in time. This rate of change is the difference between the 
rate at which new organisms are created, and the rate at which older 
organisms die. (For some populations, there might also be contributions 
from immigration and emigration.) 


To illustrate the ideas behind the way that we will model the change in a 
population size over time, let’s consider a population of bacteria in a Petri 
dish. Suppose that at a particular point in time the population size is 
100000 bacteria, and the rate of change of the population size is 30 000 
bacteria per hour. Now consider a few hours later, and suppose that the 
population size has doubled, to 200000 bacteria. Since there are now twice 
as many bacteria available to create more bacteria and to die, you would 
expect the rate of change of the population size to be twice what it was — 
that is, you would expect it to be 60000 bacteria per hour. In general, if 
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the population size is multiplied by some number, then you would expect 
its rate of change to be multiplied by the same number. 


In fact you would expect this property to hold not only for increasing 
population sizes, but also for decreasing ones. For example, consider a 
population of wild birds whose size is decreasing over time. If the 
population size decreases to half of what it was, then you would expect the 
rate of change of the population size to also decrease to half of what it 
was, since there are only half as many birds available to die and to create 
new birds. In general, if the population size is multiplied by some number 
(which can be less than 1), then you would expect its rate of change to be 
multiplied by the same number. 


If two quantities are related in such a way that whenever one of the 
quantities is multiplied by a number the other quantity is multiplied by the 
same number, then we say that the quantities are proportional (or 
directly proportional). 


So, from the discussion above, you would expect the following modelling 
assumption to hold for a population size that changes over time. 


Modelling assumption for a population size 


The rate of change of the population size is proportional to the 
population size. 


This modelling assumption does indeed work well for many populations, 
particularly over periods of time that are not too long, and for population 
sizes that are not too large. If a population size becomes too large, then 
resources such as food may start to be in short supply, so the modelling 
assumption above might not apply. 


If a quantity is proportional to another quantity, then the relationship 
between the two quantities is described by an equation of the form 


first quantity = K x (second quantity), 
where K is a constant. 


So the modelling assumption above tells you that the relationship between 
the population size and the rate of change of the population size is given 
by an equation of the form 


rate of change of the population size = K x (population size), (17) 
for some constant K. 


The constant K is called the proportionate growth rate of the 
population size. This terminology is used because, by equation (17), the 
constant K is equal to 

rate of change of the population size 


’ 


population size 
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so it is a measure of how fast the population size is increasing or 
decreasing relative to its size. For example, if the proportionate growth 
rate is K = 0.05 per year, then, by equation (17), at any point in time the 
population size is growing at the rate of 0.05 of its current size per year. In 
other words, it is growing at the rate of 5% per year. Notice that the units 
for a proportionate growth rate are ‘per year’ or ‘per hour’, or similarly for 
any other measure of time. These units can be written as year~', or 
hour~!, and so on, as appropriate. 


The larger the value of the proportionate growth rate K for a population 
size is, the more rapidly the population size will increase. A negative 
proportionate growth rate corresponds to a decreasing population size, and 
the larger its magnitude is, the more rapidly the population size will 
decrease. A proportionate growth rate of zero corresponds to a population 
size that does not change. 


Let’s write equation (17) in symbols. Let P be the population size at 
time t. Then the rate of change of P with respect to t is 


dP 
dt ` 
So equation (17) can be written as 
dP 
— = KP 


where K is the proportionate growth rate for the population. This 
equation is a differential equation that relates the time t and the 
population size P. 


We can solve this differential equation to find the relationship between t 
and P. To find the precise relationship, with no arbitrary constant, we 
need an initial condition: we need to know the population size at a 
particular point in time. Let’s assume that when t = 0 the population size 
is Po. When solving the differential equation, we can assume that P > 0, 


since P is a population size. So we have the following initial value problem: 
P 
"i = KP (P>0), where P= P) when t = 0. 


In the next activity you are asked to use the methods that you learned in 
Section 2 to solve this initial value problem. 


Activity 16 Solving the population size initial value problem 


Consider the differential equation 


dP 
i EP (P >0), 


where K is a constant. 


(a) Use separation of variables to show that the general solution of this 
differential equation is 


P = Ae*t, 
where A is an arbitrary positive constant. 
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(b) Show that the particular solution that satisfies the initial condition 
P = Py when t = 0 
is given by 
P = Pek. 


Figure 4 shows some typical graphs of the solution P = Pye*! found in 
Activity 16. The case K > 0 corresponds to population growth, K < 0 
corresponds to population decline and K = 0 gives a static population size. 


PA 


Figure 4 The graph of P = Poe*' for K > 0, K =0 and K <0 


Figure 5 illustrates that the larger the magnitude of the proportionate 
growth rate K in the equation P = Pye*', the faster the population size P 
grows or decays. 


P K=2 /,K=1 


Po 


Figure 5 The graph of P = Poe** for different values of K 


In fact, Activity 16 shows that, with the modelling assumption that you 
have seen, the variation in the size of a population over time is given by an 
exponential growth or decay function (or a constant function, in the 

case K = 0). Remember from MST124 Unit 3 that an exponential 
growth or decay function is a function of the form y = ae**, where 

a and k are non-zero constants; a positive constant k gives an exponential 
growth function, and a negative constant k gives an exponential decay 
function. You saw in MST124 Unit 3 that population sizes can often be 
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modelled by functions of this form, and the discussion above explains why 
this is; it arises from the modelling assumption discussed in this subsection. 


This model for population change is called the exponential model, or 
sometimes the continuous exponential model, for population change. 
It illustrates how an equation that relates two quantities may in fact be a 
consequence of a simpler differential equation relating the quantities. 


Here is a summary of the model. 


Continuous exponential model for population change 


If a population size changes with constant proportionate growth 
rate K, and the population size at time t = 0 is Po, then the 
population size P at time t is given by 


P= Pek". 


The next activity, which is similar to some problems in MST124 Unit 3, 
shows how the equation P = Poet can be put to use. Remember that the 
notation year~! means ‘per year’. 


Activity 17 Working with the exponential model for population 
growth 


The population size of a country is currently 5.1 million, and its 


proportionate growth rate is 0.032 year~!. 


(a) Use the equation P = Pye*! to predict what the population size will 
be after each of the following time intervals. Give your answers to two 
significant figures. 


(i) 10 years (ii) 100 years 


(b) After how long will the population size reach 50 million? Give your 
answer to the nearest year. 


A population size whose change over time is given by the exponential 
model, with a particular proportionate growth rate, has the following 
property. 
The factor by which the population size increases or decreases over a 
time interval depends only on the length of the time interval (and not 
on the time at the start of the time interval or the population size at 
the start of the time interval). 


For example, if the population size increases by a factor of 1.05 (that is, 
grows by 5%) over some time interval of length one year, say, then it will 
increase by the factor 1.05 over every time interval of length one year. 
Similarly, if the population size decreases by a factor of 0.95 (that is, 
declines by 5%) over some time interval of length one year, say, then it will 
decrease by the factor 0.95 over every time interval of length one year. You 


3 Applications of differential equations 


met this property of exponential growth and decay functions in MST124 
Unit 3. 


To see why this property holds, suppose that the population size is given 
by the equation P = Pye*', where the symbols here have the same 
meanings as earlier, and consider the time interval of length T from time t 
to time t+ T, say. The population sizes at the start and end of the time 
interval are 


Pape and per) , 


respectively. So the factor by which the population size changes over the 
time interval is 
PoeK +T) PoeKteKT 
Poekt a Poe* t 
As you can see, this value depends only on the length T of the time 
interval. (The proportionate growth rate K is a constant.) 


KT 


=e 


Example 9 Calculating a growth factor 


Suppose that a population size increases according to the exponential 
model, with proportionate growth rate 0.04 year~'. By what factor 
will the population size increase over five years? Give your answer to 
three significant figures. 


Solution 


Suppose that the initial population size, at time t = 0, is Po. The 
proportionate growth rate is K = 0.04. Hence the population size P 
at time t is given by 


P= P”, 


Consider a five-year time interval from time t to time t+ 5. The 
population sizes at the start and end of this time interval are 


Poet and Pye®-94(t+5), 
respectively. 


So the factor by which the population size will increase over the time 
interval is 
Poe O4(t+5) 0-04 0.045 


2 0.055 _ 02 = 
0.04 a 0.04t =e Se S A oan 
Poe E 


Hence the population size will increase by a factor of 1.22 (to 3 s.f.) 
over five years. 


So the population size in Example 9 will increase by about 22% over any 
time interval of five years. 
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Note that the working in Example 9 can be made a little shorter by using 
the result that the factor by which the population size changes over a time 
interval of length T is e*", where K is the proportionate growth rate. 
This formula was mentioned immediately before the example. However, to 
save having to remember this formula, you can just use the equation 

P = Pet, as was done in the solution to Example 9. 


Activity 18 Calculating a proportion remaining of a decreasing 
population size 


Suppose that a population size decreases according to the exponential 
model, with proportionate growth rate —0.03 year~!. What proportion of 
the population size will remain after a decade? Give your answer to 

two significant figures. 


You have seen that if a population size changes according to the 
exponential model, with proportionate growth rate K, then over any time 
interval of length T, the population size is multiplied by the factor eT. 


It follows that if the population size is increasing, then the time T that it 
takes to double is fixed and given by the equation 
eT _ 2. 


This time T is known as the doubling time of the population size. 


Similarly, it follows that if the population size is decreasing, then the 
time T that it takes to halve is fixed and given by the equation 
pet i 


z: 
This time T is known as the halving time or half-life of the population 
size. 


You can rearrange each of these equations to express T in terms of K, 
by first taking the natural logarithm of each side. This gives the results 
below, which you met for exponential growth and decay functions in 
general in MST124 Unit 3. 


Doubling and halving times 


Suppose that a population size changes with proportionate growth 
rate K. 


e Ifthe population size is increasing (so K > 0), then its doubling 


tj -m2 
ime is —. 
K 


e Ifthe population size is decreasing (so K < 0), then its halving 
ln 2 
time is Bees 


K 
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Activity 19 Finding doubling and halving times 


Find, to three significant figures, the doubling time for the population 
described in Example 9, and the halving time for the population described 
in Activity 18. 


3.2 Radioactive decay 


Radioactive substances are composed of atoms that ‘decay’ spontaneously; 
that is, they are transformed into other types of atoms, releasing energy 
and subatomic particles (both often called radiation) in the process. The 
rate at which this occurs varies according to the radioactive substance 
concerned. 


For example, the radioactive element ruthenium-106 slowly decays into 
another element, rhodium-106. Each atom of ruthenium-106 has a known 
probability of decay in any given length of time, which is a constant 
throughout time. In fact each atom of ruthenium-106 has almost exactly a 
50% chance of decay in a year. So if we have 10000 atoms of 
ruthenium-106, about 5000 will decay in a year, and if we have 100000, 
about 50000 will decay in a year. 


Because the number of atoms present in even one gram of ruthenium-106 
is enormous (about 5.7 x 10?! atoms), it is sensible to use a continuous 
model to describe this situation, measuring the amount in grams, rather 
than in numbers of atoms. 


With such a model, we can consider the rate of change of the amount of 
ruthenium-106 at any particular moment in time. This rate of change will 
be proportional to the amount of ruthenium-106 present. This is because, 
for example, if the amount of ruthenium-106 is halved, then there are only 
half as many atoms to decay (each with the same probability in any given 
length of time), so the rate of change of the amount of ruthenium-106 will 
also halve. 


Other radioactive atoms have different probabilities of decay in a given 
time interval, but these are also a constant throughout time and hence a 
similar model will hold. In general, the rate of decay of any radioactive 
substance can be modelled by the following relationship. 


Radioactive decay 


The rate of change of the amount of a radioactive substance is 
proportional to the amount of the substance present. 
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So we have the following equation, where K is a constant: 
rate of change of amount of 
( the amount of the = Kx ( the aati | . 
radioactive substance substance 
Since the amount of the radioactive substance is decreasing, the 
constant K in the equation above is negative, and we usually write it 
as —k, where k is a positive constant, called the decay constant. Like a 


proportionate growth rate, a decay constant is measured in units such as 
year !, hour~', and so on. 


Let’s write the equation above in symbols. Let m be the amount (mass) of 
the radioactive substance at time t. Then the rate of change of the amount 
of the substance at time t is dm/dt. 
So the equation can be written as 
dm 
dt 
where k is the decay constant. 


= —km, 


To find a particular solution of this equation we require an initial 
condition. If the amount of the radioactive substance present at time t = 0 
is mo, then the initial condition is 


m = mo when t = 0. 


So we have the following initial value problem: 


dm 


= Thm (m>0), where m = mo when t = 0. 


The condition m > 0 holds because the amount of the substance is always 
positive. 


This initial value problem has the same mathematical form as the initial 
value problem for the population change model (with —k replacing K). 
Hence its solution can be written down immediately, by comparison with 
the exponential model for population change. The solution is therefore 
m= moe". 


So we have the following result. 


Exponential model for radioactive decay 


If a radioactive substance decays with decay constant k, and the 
original amount of the substance is mo, then the amount m of the 
substance at time t is given by 


im = nem 


The graph of the solution m = mge~*, shown in Figure 6, illustrates that 
this solution has the characteristics we would expect of radioactive decay. 
The amount of the radioactive substance decreases, initially quite quickly 
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and then more slowly, until eventually there is practically none of the 
original substance left. 


MLA 


Mo 


m= moe Ft 


œY 


Figure 6 The graph of m = moe™** for k > 0 


Different radioactive substances have different values for the decay 
constant k. The larger the value of the decay constant k, the faster the 
substance decays. This is illustrated in Figure 7. 


Figure 7 The graph of m = moe™** for several values of k 


Because the solution m = moe™** is an exponential decay function, the 
properties of exponential decay functions discussed in the last subsection 
hold for radioactive decay. In particular, every radioactive substance has a 
half-life, which is the time that it takes for the amount of the substance 
to halve, and the relationship between the decay constant and the half-life 
is as follows. 


Half-life of a radioactive substance 


If a radioactive substance decays with decay constant k, then its 
half-life T is given by 
In2 


T === 
k 


Although the speed of decay of a radioactive substance can be described in 
terms of the decay constant k, it is more usual to use the half-life T. 
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Ernest Rutherford 
(1871-1937) 


Skull of a smilodon 
(sabre-toothed cat). Carbon 
dating shows that this exotic 
cat, larger than a modern 
tiger, died about 10000 years 


ago. 
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Different radioactive substances have a very wide range of half-lives, from 
billions of years down to fractions of a second. 


Activity 20 Using a half-life 


Silicon-31 has a half-life of 2.62 hours. 


(a) Find the value of the decay constant k for silicon-31, to three 
significant figures. 


(b) Hence calculate, to two significant figures, the proportion of a sample 
of silicon-31 that will still be present after 6 hours. 


Ernest Rutherford was a New Zealand-born British physicist of great 
distinction. He discovered the concept of radioactive half-life and was 
the first to show that radioactivity was caused by the decay of heavy 
atoms into lighter atoms, for which he was awarded the 1908 Nobel 
Prize in Chemistry. 


Carbon dating 


One useful application of radioactivity is for dating archaeological finds of 
a biological origin. Carbon-14 is a radioactive element present in the 
atmosphere that decays (to nitrogen) with a half-life of about 5570 years. 
This decay is balanced by the production of new carbon-14 in the 
atmosphere by cosmic rays, so that the proportion of carbon-14 in the 
atmosphere is approximately constant. It is absorbed by living organisms 
when they are alive, but is no longer replenished after the organism dies, 
so the amount of carbon-14 in the remains of the organism declines as the 
carbon-14 decays. By measuring the proportion of carbon-14 in a dead 
organism and comparing it to the proportion in a living one, the number of 
years since death can be estimated. This technique is known as carbon 
dating or radiocarbon dating. 


Activity 21 Using carbon dating 


Suppose that the bones of an animal are found in an archaeological dig 
and analysis produces the estimate that 85% of the original amount of 
carbon-14 is still present in the bones. Taking the half-life of carbon-14 to 
be 5570 years, find the approximate age of the bones to the nearest 

100 years. 
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3.3 Newton’s law of cooling 


When you make a hot cup of tea and leave it on your desk, the tea cools 
down, eventually reaching the ambient temperature of the room. 
Furthermore the colder the room is, the faster the tea cools down. This 
observation is the foundation of the following law. 


Newton’s law of cooling 


The rate at which a body cools is proportional to the temperature 
difference between the body and the surrounding medium. 


A thermal image of a recently 
poured cup of tea. Orange 
indicates the warmest areas, 
blue and purple the coolest. 


In our example the body is the cup of tea and the surrounding medium is 
the air in the room. 


Empirical observation has shown that Newton’s law of cooling holds, to a 
good approximation, in a wide variety of situations where the body and 
the surrounding medium are both good conductors of heat. 


Let’s write Newton’s law of cooling in symbols. Let T be the temperature 
of the body (such as the cup of tea) at time t, and let Ta denote the 
ambient temperature of the surrounding medium (such as the air), where 
the subscript ‘a’ is short for ‘ambient’. Then the rate at which the body 
cools is dT/dt, and the temperature difference between the body and the 
surrounding medium is T — Th. 


Hence, according to Newton’s law of cooling, we have 


dT 
rT K(T —T,), 


where K is a constant. Since the temperature of the body is decreasing, 
the value of dT/dt is always negative. Also, the ambient temperature is 
always lower than the temperature of the body, so T — Ta is always 
positive. Therefore K is a negative constant, and hence we can write 

K = —X, where À is a positive constant. The equation above then becomes 


dT 

n 
The positive constant À is dependent on both the body and the 
surrounding medium. 


-A(T -T,). 


If we assume that the surrounding medium is large in comparison to the 
hot body (as in the case of the air surrounding the cup of tea), then it will 
hardly heat up as the body is cooling down, so we can assume that Th is a 
constant. Under these assumptions we can solve the differential equation 
above using separation of variables. In doing so we can assume that 

T — Ta > 0, a fact mentioned above. 


To find a particular solution of the differential equation, we need an initial 
condition. If the initial temperature of the hot body is Tọ, then the initial 
condition is 


T = To when t = 0. 
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In the next activity you are asked to solve the differential equation that 
arises from Newton’s law of cooling, with the initial condition above. 


Activity 22 Solving Newton's law of cooling 


Consider the differential equation 
dT 
a 

which arises from Newton’s law of cooling. 


=A — Ta) (T-T, > 0), 


(a) Use separation of variables to show that the general solution of this 
differential equation is 


T=T,+Ae™, 
where A is an arbitrary positive constant. 
(b) Show that the particular solution that satisfies the initial condition 
T = To when t= 0 


T =T; + (To — Taje ™™. 


(c) Suppose that a body cools according to Newton’s law of cooling. Use 
the particular solution found in part (b) to find the eventual 
temperature of the body after a long time. 


Here is a summary of the solution determined in Activity 22. 


Solution of Newton’s law of cooling 


If a body cools in a surrounding medium with constant 
temperature Ta, and the temperature of the body at time t = 0 is Tọ, 
then the temperature T of the body at time t is given by 


T E E lene 


where À is a positive constant. 


Activity 23 Determining a time of death 


The temperature of a healthy human being is 37°C. When a person dies, 
the body starts to cool, and its rate of cooling depends on the surrounding 
ambient air temperature. 


Suppose that a dead body is found and a pathologist is called in to 
determine the time of death. The pathologist finds that the room in which 
the body is located is kept at a constant temperature of 20°C. 
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She takes the temperature of the body when she first arrives and finds it to 
be 30°C. She takes another temperature reading an hour later and finds it 
to be 28.5°C. 


(a) By finding the values of the constants in the solution of Newton’s law 
of cooling that apply to this situation, show that the temperature T of 
the body as a function of time t, taking t = 0 to be the time at which 
the pathologist arrived, is 


T = 20+ 10 x 0.85". 


(b) How long before the pathologist arrived did the person die? Give your 
answer to the nearest quarter of an hour. 


4 Solving linear differential equations 


This section presents a third method of finding solutions to first-order 
differential equations. The method, called the integrating factor method, 
applies only to a particular form of equation known as a linear differential 
equation. 


4.1 Linear differential equations 


A first-order differential equation involving the independent variable x and 
the dependent variable y is linear if it can be expressed in the form 


d 

T +g(2)y = h(a), (19) 
x 

where g(x) and h(x) are given functions. 


The linear first-order differential equation (19) is homogeneous if 
h(x) = 0 for all z, and inhomogeneous or non-homogeneous otherwise. 


The term linear here refers to the dependent variable y. Although any 
functions g(x) and h(x) of the independent variable x are allowed in the 
differential equation, the only occurrences of y allowed are those in the 
terms dy/dx and g(x)y. This means that terms involving y?, y’, siny, and 
so on, are not allowed. 


So, for example, the differential equation 
dy 


2 3 
——ry=2 
dx u 
is linear, with g(x) = —x? and h(x) = x, whereas the differential equation 
dy 2 
— == 
dx 


is not linear, due to the presence of the term y?. 


A human hand heating the 
surrounding still air 
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Activity 24 Identifying linear differential equations 


State which of the following differential equations are linear, identifying 
the functions g(x) and h(x) in each case that is linear. 
dy 3 5 d dz 


(a) +ry=r (b) q Tesina (c) T 3z 


d d 
(d) Y jyt (e) r y= (£) (1422) + ory = 32? 
dx dx 


Note that if g(x) = 0, then equation (19) reduces to 


st = h(a), 
which is the directly integrable type considered in Subsection 1.2. So, 
directly integrable differential equations are also linear differential 
equations. Since we have already considered this case in some detail, for 
the remainder of this section we will concentrate on the case where 


g(x) #0. 


4.2 Constant-coefficient linear differential 
equations 


One of the simplest examples of a linear differential equation is the case 
where h(x) = 0 (so the equation is homogeneous) and g(x) = A, where A is 
a constant. In this case, equation (19) becomes 


dy 
— + Ay =0. 
dx Tey 
We dealt with this differential equation in some detail in the previous 


section (in the equivalent form dP/dt = KP). It has solution y = yoe~4* 
where yo is a constant. 


’ 


Now consider the case where h(x) 4 0 (that is, the differential equation is 
inhomogeneous) and g(x) = A, where A is a constant. This gives the form 


d 
IY 4 Ay = h(a), (20) 
dx 
which is called a constant-coefficient first-order linear differential 
equation. This is because in the context of the linear differential 


equation (19), the function g(x) is often called the coefficient of y. 
We can use a trick to solve differential equations of this form. Here’s the 
idea. First multiply both sides by e^”: 
d 
etn + AeA? y = eA h(a). (21) 
z 


Now, by the product rule for differentiation, 
d 


x zd x 
a y) = e^ = + AeAy, 


4 Solving linear differential equations 


So the left-hand side of equation (21) is equal to the derivative 
d 


qe) 

and hence equation (21) can be written as 
d x T 
Ta y) = e4*h(a). 


This equation can be solved by integrating both sides with respect to x, 
which gives 


eAty = fne) dz. 


Hence the general solution of differential equation (20) is 


y=e 4 ( f eA h(a) ar) . (22) 


So, if we can find the integral on the right-hand side, then we can solve 
differential equation (20). Here is an example where the integral can be 
found by using integration by parts. 


Example 10 Solving a simple linear differential equation 


(a) Solve the differential equation 


(b) Check that your solution satisfies the differential equation. 
Solution 


(a) ®. The differential equation has the form of equation (20), with 
A= 2 and h(x) = x. So the solution is given by equation (22). @ 


By equation (22), the general solution is 


y=e (J re??? ar) : 


@. The integrand is a product of two expressions, x and e°”. You 
can integrate the second expression, e°”, and differentiating the 
first expression, x, makes it simpler. So try integration by 

parts. @ 


Integrating by parts gives 


y= e722 (gee = fe ar) 
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®. Do the integration inside the brackets, remembering to 
include the constant of integration inside the brackets. © 


So, the general solution is 
y= ze" (ze = se" +c) 
=} (r-i +e”), 
where c is an arbitrary constant. 

(b) Differentiating the general solution we get 
dy 
dz 

Substituting this and the general solution into the left-hand side 
of the differential equation gives 


= + (1 — Jee =) : 


— + 2y = } (1 — 2ce”) +2 (x — $+ ce”) 
i- ce” tr- itce” 
= I, 


which is equal to the right-hand side. So y = 3 (x — 2 + ome) 
satisfies the differential equation. 


Here is a similar example for you to try. 


Activity 25 Solving a simple linear differential equation 


(a) Find the general solution of the differential equation 


dy 2a 
de 7 28 : 


(b) Check that your solution satisfies the differential equation. 


In the next subsection you'll see that equation (22) can be adapted to 
solve the case where the coefficient g(x) is not constant. 


4.3 The integrating factor method 


The key to solving the constant-coefficient linear differential equation was 
to multiply both sides by e4*. This allowed us to integrate the left-hand 
side of the equation. For this reason e^” is called an integrating factor 
for the differential equation. We now generalise this trick to solve the 
general linear first-order differential equation (19). 


4 Solving linear differential equations 


Integrating factor method 
To find the general solution of a differential equation of the form 


d 

E t+ 9(a)y = h(a), (23) 
determine the functions g(x) and h(a), and then perform the 
following steps. 


1. Write down the integrating factor 


plz) =exp( f g(a)ar). (24) 


Find the integral, if possible, but don’t include a constant of 
integration. 


2. Write down the general solution, given by the formula 


me => ( / BAG) ar) , (25) 


Find the integral, if possible, and include the constant of 
integration inside the brackets. 


Remember that exp(x) is just another way to write e”, so equation (24) 


can be written as 
p(z) = efs, 


You will see why the integrating factor method works later in the section. 
Notice, however, that it is a generalisation of the method that was used for 
the constant-coefficient linear case in Subsection 4.2. In that case 

g(x) = A, so the integrating factor given by equation (24) is 


=o ( f Ada) e 


This is exactly the integrating factor that was used in Subsection 4.2. The 
general solution given by equation (25) is therefore 


y=e 4 ( f eA? h(x) ar) 


exactly the same as equation (22). 


As with the separation of variables method, when you use the integrating 
factor method it is not always possible to find the necessary integrals. 
However, this is always possible in the examples and activities in this unit. 
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— Use the integrating factor method to find the general solution of each 
of the following differential equations. 


dy oe dy il 
(a) qe Ye sing (b) ae oe (x > 0) 
Solution 


(a) &. Verify that the given differential equation is linear. ® 
The given differential equation has the form of equation (23) with 
e= =l wand Mo =e" cme. 
so we use the integrating factor method. 


@. First calculate an integrating factor, without including a 
constant of integration. © 


By equation (24), the integrating factor is 
p(x) = exp (fn ar) = exp(—x) = e”. 


®. Calculate the general solution, including the constant of 
integration inside the brackets. ©& 


The general solution given by equation (25) is 


y= = (S e sima) ar) a ( [ sina) l 


Doing the integration, we obtain 
y = e” (— cosg + c), 
where c is an arbitrary constant. This is the general solution. 
(b) ®. Verify that the given differential equation is linear. ® 


We rearrange the differential equation so it is in the form of 
equation (23): 


aly i il 

W Pk 26 

da i > 
This now has the form of equation (23), with 


1 il 
g(a) =-— and h(x) =-=, 


so we use the integrating factor method. 
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@. First calculate the integrating factor, without including a constant 
of integration. ® 


By equation (24), the integrating factor is 


= exp (— nal) 
= exp (— ln z) (since x > 0) 
= exp (In (a=) 

1 


@. Calculate the general solution, including the constant of 
integration inside the brackets. © 


The general solution, given by equation (25), is 


EU OCS) 


Il II 
8 8 
ee 
Sli | 
SF = 
SEPE 
eee, b| = 
a 
8 
Sense 


= 1+ ee. 


where c is an arbitrary constant. This is the general solution. 


Note that the differential equation in part (b) of Example 11 can also be 
solved by the method of separation of variables, if you assume that y Æ 1 
(see ‘A note of caution’ on page 176). You might like to check that you get 
the same answer if you solve it this way (except that one particular 
solution is omitted because of the condition y Æ 1). 


Here are some differential equations that you can solve by the integrating 
factor method. 


Activity 26 Using the integrating factor method to find general 
solutions 


Find the general solution of each of the following differential equations. 
dy dy 


(a) so TtT (b) ve sya (x > 0) 
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Activity 27 Using the integrating factor method to find particular 
solutions 


Solve each of the following initial value problems. 
du 
(a) a 


xu, where u(0) =2 


(b) 2 +2y=t (t>0), where y(1)=1 


Activity 28 Newton's law of cooling: change of ambient 
temperature 


In Subsection 3.3 you met Newton’s law of cooling, given by the 
differential equation 

dT 

dt 
where A is a positive constant. 


The solution of this differential equation describes the temperature T of a 
body at time t, when it is placed in a surrounding medium of ambient 
temperature Ta. In that subsection we assumed that the ambient 
temperature Tą was a constant. However, if Tẹ, varies reasonably slowly 
with time, the equation remains valid and it can be rearranged to give 

Z +AT = AN, (t), (27) 


which is a first-order linear differential equation. 


(a) Suppose that the ambient temperature Ta is a function of time that 
decays exponentially with time; that is, 


T(t) = Tao e™, 


where k and Thao are positive constants, and 0 < k < à. This equation 
models a surrounding medium that has ambient temperature Tao at 
time t = 0, and then cools to temperature 0 after a very long time. 
Then equation (27) becomes 


dT 
— AT = Aloe". 
q tÀ À Oe 


Show that the general solution of this equation is 


where A is an arbitrary constant. 


(b) Check that when k = 0 the function T,(t) is a constant and the general 
solution reduces to that obtained in Activity 22. 


4 Solving linear differential equations 


As always, when you have solved a differential equation it is a good idea to 
check that the solution you obtain satisfies the original differential 
equation. As before, you can do this by substituting your solution and its 
derivative into the differential equation. 


Example 12 Verifying a solution 

Verify that the solution to Example 11(a) satisfies the differential 
equation in that example. 

Solution 


The differential equation was 
d 
ae Y= e Sint 


daz 
and its general solution was calculated to be 


y = e” (— cosx + c). 
Differentiating this general solution using the product rule gives 


W aa oT 
n gz cose te) + g Eo cos xz + c) 


= e” sing + e” (— cosg + c). 
Substituting for y and dy/dz in the left-hand side of the differential 
equation gives 
dy 


La e” sin x + e”(— cosx + c) — e” (— cosx + c) 
a 


= E sim ae. 


which is equal to the right-hand side, as required. 


Here’s an example for you to try. 


Activity 29 Verifying another solution 


Verify that the solution to Example 11(b) satisfies the differential equation 
in that example. 


4.4 Why the integrating factor method works 


If you would like to know why the integrating factor method works, then 
you can read the following justification. Otherwise, skip this material and 
go on to the next section. 
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The general first-order linear differential equation is of the form 


d 
T + 9(e)y = hla), (28) 
x 
and we show that its general solution is 
1 
= — x)h(x) dz }, 29 
y= ([ (wine) ar) (29) 
where 


ple) = exp ( f (a)ar). 


Here p(x) is the integrating factor. 


The integrating factor p(x) is chosen to be of this form because, when it is 
differentiated with respect to x, we obtain the following equation: 


dp 
ac 30 
P = ple)gla) (30) 
To see why this equation holds, we set p(x) = e”, where u = f g(x) dz, and 
then use the chain rule for differentiation to give 

dp dg, dt y B 

d ue Jz = e"g(x) = p(x) g(x), 
as required. 


By the product rule and equation (30), 


= paj) = ola) + Py 
= ple) Z + paga). (31) 


Now, multiplying both sides of equation (28) by p(x) gives 


pa) E + ple)sle)y = pe)hle). 


The left-hand side of this equation is the same as the right-hand side of 
equation (31), so 


It follows that 
(oly = | poha) ax. 
Dividing both sides by p(x) gives the general solution as 


y= ( [roras], 


which is equation (29), as required. 


5 Summary and consolidation 


Finally, you may have been wondering whether we would obtain further 
solutions to this differential equation by including a constant of integration 
when we find the integrating factor. But if we did this, then the 
integrating factor would be of the form 


ne een (fate) ice c) sdi ¢ a ar) , 


where C is a constant, so the new integrating factor would simply be a 
constant multiple of the original one, and in equation (29) this constant 
multiple would be cancelled out. 


5 Summary and consolidation 


It’s time to summarise and review the main results of this unit so far. You 
have seen three types of first-order differential equation: 


dy _ dy 


an f(z), == = f(x)g(y), LF g(x)y = h(x). 


The first type is directly integrable and is solved by direct integration, the 
second type is separable and is solved by separation of variables, and the 
third type is linear and is solved by the integrating factor method. Also, 
remember that the first type is a special case of the other two. 


You need to be able to recognise each type of differential equation and 
know its method of solution. It is important to realise that often variables 
other than x and y are used, and that a differential equation may need to 
be rearranged in order to put it into one of these forms. 


Here is an activity to check your ability to recognise each of these three 
types of differential equations, followed by several activities that are 
intended to strengthen your ability to solve them. 


Activity 30 /dentifying differential equations 


Which method would you use to try to solve each of the following 
differential equations? If more than one method can be used, you should 
list all the methods. (You are not asked to solve these equations.) 


(a) OY og cB ca (b) C T 


dg dg 

dv dy 

2 = 1 2) Z 4 2xy = 1 — x? 
(e) ga (d) ( tel + TY z 

dy 2y du 2 2 

Å= M = =< I — — 4 
(e) eA y+t (£) Jz te + 2(u* — x) 
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Activity 31 An application: a melting snowball 


The rate of change of the radius, r, of a melting snowball with respect to 
time, t, can be modelled by the differential equation 

dr 

dt 
where k is a positive constant. It is assumed that the snowball is spherical 
at all times and that the temperature of the air around it is a constant 
(above 0°C). The radius r of the snowball is measured in centimetres and 
the time t is measured in days since melting started. 


=—k (r 2), 


(a) Find the general solution of this differential equation. 


(b) The snowball originally has a radius of 100cm. Find the particular 
solution that describes its radius as a function of time. 


(c) The snowball is found to have lost half its initial volume after 2 days. 
Show that k = 50 (1 — ont), (Remember that the volume V of a 
sphere is given in terms of its radius r by V = nr.) 


(d) How many days does it take for the snowball to disappear completely? 
Give your answer to the nearest day. 


Activity 32 Finding a general solution and a particular solution 
(a) Find the general solution of the differential equation 
dz 
— =tcost. 


dt 


(b) Find the particular solution that satisfies the initial condition z = 2 
when t = 0. 


Activity 33 Solving an initial value problem 
(a) Solve, in explicit form, the initial value problem 


d 
= =e (x >0), where xz = 2 when t= 0. 


6 Other methods of solving differential equations 


Activity 34 Solving an initial value problem and checking the 
solution 
(a) Solve, in implicit form, the initial value problem 
dz t 
dt cosr+a2 ei), Nase aly) 


(b) Verify that the particular solution that you have found satisfies the 
differential equation. 


Activity 35 Solving an initial value problem 


Solve the initial value problem 


d 


7 1—e%y, where y(0) = 3. 


Activity 30 Finding a general solution and a particular solution 


(a) Find the general solution of the differential equation 


dy Qxry 
=f f 
dz z2? +1 


(b) Find the particular solution of the differential equation that satisfies 
y(0)=3. 


6 Other methods of solving differential 
equations 


This final section describes various other approaches to finding solutions of 
first-order differential equations and to understanding the behaviour of 
these solutions. 


Consider the general first-order differential equation of the form 


dy _ 
dz — f(x,y), 


where f(x,y) is some given function of x and y. As you have seen, for 
certain special types of function f(x,y), you can solve the differential 
equation, obtaining a general solution for y in either explicit or implicit 
form, and then obtain particular solutions by using initial conditions. 
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If such a solution to a differential equation can be obtained, then the 
differential equation is said be solvable analytically and the solution is 
said to be an analytic solution. 


There are more elaborate techniques, which are beyond the scope of this 
module, for finding analytic solutions when f(x,y) has certain more 
complicated forms. Indeed, the module computer algebra system can find 
analytic solutions to a much wider variety of first-order differential 
equations than you have seen here. Sometimes such solutions are expressed 
in terms of functions that are familiar, like sin, cos, exp, In, and so on, and 
in other cases the solutions are in terms of more obscure functions. 


However, there are many functions f(x,y) for which the differential 
equation dy/dx = f(x,y) has no analytic solutions. For example, the 
differential equation 


= = sin(xy), 


dx 


has no analytic solutions in explicit or implicit form. Such differential 
equations occur frequently in the mathematical sciences, and are probably 
more common than differential equations with analytic solutions! 


Mathematicians and scientists have developed a number of techniques to 
explore the solutions of such differential equations. These include: 


e finding approximate analytic solutions 
e plotting direction fields 
e calculating numerical solutions. 


The first of these techniques consists of a wide range of methods for 
constructing approximate analytic solutions to differential equations, which 
are close to exact analytic solutions. This is a vast area of applied 
mathematics, well beyond the scope of this module, so we will say no more 
about it. 


The second technique, plotting direction fields, is described in 

Subsection 6.1. The third technique, calculating numerical solutions, is 
discussed briefly in Subsection 6.2. Then, in Subsection 6.3, the use of the 
computer for solving differential equations is covered. 


6.1 Direction fields 


A direction field is a graphical representation of a differential equation. 
To appreciate the idea, consider the simple differential equation 


dy > 
dr’ 


from Subsection 1.2, which has general solution 
y = 42° +c, 


where c is an arbitrary constant. This general solution consists of an 
infinite family of functions, some of whose graphs are shown in Figure 8. 


6 Other methods of solving differential equations 


C= 43 
C= il 

C= 2/3 
e= iis 


2j 0 lee 


Figure 8 The graphs of y = ir’ +c, for c = Z, 2, 1 and $ 
At each point (x,y) on one of these graphs, the equation dy/dx = x? must 
hold; that is, the gradient of the graph is z?. For example, at the point 


(x,y) = (1,1) the gradient is z? = 1? = 1. 


This property is also true for a general differential equation of the form 
dy/dx = f(x,y). That is, each solution of the differential equation is a 
function, y = y(x), with the property that at each point (x,y) on the 
graph of the function the gradient is f(x,y). When these gradients are 
depicted graphically at many points of the plane, we obtain what is called 
a direction field. 


This is done for the simple case of f(x,y) = x? in Figure 9. Here a grid of 
points has been chosen and at each point (x,y) of the grid a short line 
segment has been drawn centred on (x,y), with the gradient x”. The idea 
is that if the spacing between the grid points is not too large, then a good 
picture is obtained of the shape of the graphs of the solutions to the 
differential equation. 


In this figure, the graph of one of the particular solutions, y = ir’ + 3, has 


also been plotted. As you can see, the short line segments near this graph 
are close to being tangent to it. 


NE Sao N ee Nar Ss 


NS NEON NNN NES 


SS SSS SSS NS 
VONA Nee Na NG NN 
NON NY WN XN NNN 


SVN So oa i Ne ON 


Figure 9 A direction field for dy/dx = x, together with the particular 
solution y = 42° + 3 
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The distribution over the plane of such short line segments with the 
appropriate gradients is called a direction field for the differential 
equation. Even when you can’t find the solution to a differential equation 
dy/dx = f(x,y) directly, you can always plot a direction field for the 
differential equation. 


Figure 10 shows a direction field for the differential equation 


dx 


which can’t be solved analytically. You can see that the gradients of the 
line segments plotted do correspond to the numerical values of 

f(x,y) = sin(xy). For example, these gradients are 0 when x = 0 (that is, 
along the y-axis) and also when y = 0 (that is, along the x-axis), and the 
gradient at the point (x,y) = (1,1) is sin 1 =0.841.... 


= sin(xy), 


AN NN, yt ae ee a a 
- ȚS N S ~n - - 24 Oe = SON SS ~ 
ee a NON 2 a a A a ee N ON ON 
SON O NES L a a N A CA a a ea ENS 
W S SS. a A 7 ais 
M ‘N SS SS ES Pati te 7 J rn E a a 
A N nN NO ee 7 7 7 7 1 
Se = a ee, 7 7 
~ [x ~ ~ = =| - - =| C - Pa a 1 a 
eee iia e ë e 
i i i > a 


Figure 10 A direction field for dy/dx = sin(xy), and the particular 
solution satisfying the initial condition y(0) = 0.5 


The graph of any particular solution to the differential equation must have 
a gradient at each point that matches the gradient provided by the 
direction field. You can see that the gradients at points on the graph of 
the particular solution plotted in Figure 10 match the gradients of the 
nearby line segments very closely. Thus the direction field gives a good 
intuitive impression of how solutions of this differential equation behave. 
If the short line segments were plotted more closely together, then this 
impression would be even clearer. 


So, a direction field is a graphical representation of a differential equation, 

which gives you an idea of what the solutions are like. Direction fields are 

most often used for differential equations that can’t be solved analytically, 

but they are also useful even if an analytic solution can be obtained, giving 
a qualitative understanding of the behaviour of the general solution. Such 

graphical representations are used extensively in more advanced studies of 
differential equations. 


In principle, direction fields can be plotted by hand, although it is a 
tedious task. In Activity 37, you will see how the module computer algebra 
system can be used to plot direction fields. 


Learning outcomes 


6.2 Numerical solutions of differential equations 


You have seen that the direction fields introduced in Subsection 6.1 give 
useful qualitative information about the solutions of a differential equation. 
This geometric approach can be developed further in order to use a 
computer to produce quantitative information about the solutions of the 
differential equation, in the form of numerical solutions. 


A numerical solution of the differential equation dy/dzx = f(x,y) 
consists of a sequence of points in the x, y-plane. These points are chosen 
in such a way that taken together the points form an approximation to the 
graph of a particular solution of the differential equation. Roughly 
speaking, they are obtained by ‘joining up’ short line segments in a 
direction field for the differential equation. There are various ways of 
calculating such a sequence of points, the details of which are beyond the 
scope of this module, but the idea is as follows. 


We start with a given value of x and the corresponding value of y, given by 
the initial condition. The value of x is then increased by a small amount to 
a new value and an estimate is obtained for the corresponding new value 
of y by using the expression f(x,y) for the derivative dy/dx. This step is 
then repeated many times until an approximate numerical solution has 
been obtained for the required range of values of x. The more points there 
are in the sequence, and hence the closer together adjacent points are, the 
more accurate the numerical solution is likely to be. 


6.3 Solving differential equations using a computer 


In the following activity, you can learn how to use the computer algebra 
system to find analytic and numerical solutions of differential equations 
and to investigate direction fields. 


Activity 37 Solving differential equations with a computer ne 
To complete your study of Unit 8, work through Section 7 of the Computer E 
algebra guide, where you will learn how to obtain analytic and numerical 
solutions of differential equations, and how to plot direction fields. 
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Learning outcomes 


After studying this unit, you should be able to: 


recognise differential equations that can be solved by direct 
integration, and use direct integration to solve them 


recognise separable differential equations and use the method of 
separation of variables to find the general solution in explicit or 
implicit form 


recognise first-order linear differential equations and use the 
integrating factor method to solve them 


use an initial condition to find a particular solution of a differential 
equation, given the general solution 


check a solution of a differential equation in explicit or implicit form 
by substituting it into the differential equation 


understand the terms exponential growth/decay, proportionate growth 
rate and doubling/halving time, when applied to population models 


understand the terms exponential decay, decay constant and half-life, 
when applied to radioactivity 


solve problems involving exponential growth and decay 


apply differential equations and their solutions to simple models of 
real-world problems 


understand the ideas of direction fields and numerical solutions, and 
find these using the computer algebra system. 


Solutions to activities 


Solution to Activity 1 


Directly integrable differential equations have the 
form 
d(dependent variable) 


d(independent variable) 
Equation (a) is of this form. Equations (b) and (c) 
are not. 


Solution to Activity 2 


(a) Integrating both sides of the differential 
equation with respect to x gives 


y= f (2r+3)de 


=x? + 3r +c, 
where c is an arbitrary constant. 
(b) With y = x? + 3x + the left-hand side of the 


differential equation gives 


d d 
T EF (£? + 3x +c) = 22+ 3, 
which is the same as the right-hand side, so the 
differential equation is satisfied. 


Solution to Activity 3 
(a) The general solution of da/dt = sin t is 


y= fsintat 


= — cost + cC, 
where c is an arbitrary constant. 
(Check: differentiating x = — cos t + c gives 
dz/dt = sin t, as required.) 


(b) Only (iii) is a particular solution of the 
differential equation. It corresponds to c = 1. 


Solution to Activity 4 

Substitute the values x = 1 and y = 7 into the 

equation y = 5,/z +c. This gives 
7=5Vl+c=5+te, 

so c = 2. The required particular solution is 

therefore 


y = 5yz +2. 


= f (independent variable). 


Solutions to activities 


Solution to Activity 5 
(a) The general solution of dy/da = yT = x'/? is 


y= | oar 


= 2g/ 2 +e, 
where c is an arbitrary constant. 
The initial condition is y = 5 when z = 4. 


Substituting x = 4 and y = 5 into the general 
solution gives 
5=2x 7% caBie, 


soc= —. Hence the particular solution that 


satisfies the initial condition is 
y= 223/ Bis £. 


(b) The general solution of du/dxz = cos(2x) is 
y= [ costae) dx 


= į sin(2x) + c, 


where c is an arbitrary constant. Substituting 


u = —2 and x = $r into the general solution 
gives 

—2=4sin($7)+c=$+t¢, 
soc= —3. Hence the particular solution that 


satisfies the initial condition is 


u = į sin(2x) — 3. 


Solution to Activity 6 
The general solution of dz/dt = k/t (t > 0) is 


k 
v= Ea 
t 


= kln|t|+c, 
where c is an arbitrary constant. However, we have 
the restriction t > 0, so the general solution is 
x= klnt+c. 
Using the initial condition, x = 3 when t = 1, we 
obtain 
3=klnl+c=c. 


Hence c = 3, and the solution of the initial value 
problem is 


x= klnat+3. 
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Solution to Activity 7 
(a) The general solution of dm/dt = —8t is 


m= f (—8t) dt 


= —4 +c (t>0), 
where c is an arbitrary constant. (This general 
solution can be checked by differentiation.) 


(b) The initial condition is m = 200 when t = 0. 
Substituting t = 0 and m = 200 into the general 
solution gives 


200 = 0 + c, 


so c = 200. Hence the required particular 
solution is 
m = 200 — 4t? (t>0). 

(c) When m = 100, the corresponding value of t is 

given by 

100 = 200 — 4#?; thatis, t= 25. 
Hence t = 5 when m = 100, so the rocket’s fuel 
is exhausted after 5 seconds. 


Solution to Activity 8 
Equations (a), (b), (c) and (e) are of the 
appropriate form, but equations (d) and (f) are not. 


(a) f(z) =2, g(y)= cosy 


(b) fi) =teost, g(z)=1 
(c) f(t)=1, g(p) = pcosp 
(e) fy=1+y, glv)=v 


These are the most likely choices for f and g, but 
they are not the only possible answers. For 
equation (a), for example, you could have chosen 
f(x) = 2x and g(y) = $ cosy. 

(Equation (b) can also be solved using direct 
integration; this is the most straightforward method 
of solution.) 


Solution to Activity 9 


(a) The differential equation dy/dz = —x/y has a 
right-hand side of the form f(x)g(y), where 
f(x) = —zx and g(y) = 1/y, so it is separable. 
By separating the variables, we obtain 


[uav= fiaa 
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Carrying out the two integrations gives 
sy" = $x" PG; 
where c is an arbitrary constant. Since the 


left-hand side is positive when y > 0, the 
constant c must be positive. 


This implicit form of the general solution can 
be manipulated to make y the subject, giving 


Yr=xuv 2c — z2. 


However, we have a restriction of y > 0 in the 
question, so 


y = V 2c — 2x2, 
or, equivalently, 
YEN b— a", 


where 6 = 2c is another arbitrary positive 
constant. 


(The implicit form of the general solution can 
be written in the form x? + y? = b, which is the 
equation of a circle. The condition y > 0 led to 
the explicit solution y = vb — x?, whose graph 
is the upper half of this circle.) 

(b) The initial condition is y = 1 when x = 0. 
Substituting x = 0 and y = 1 into the general 
solution gives 


1=Vb-0 = vh, 


so b = 1. Hence the particular solution that 
satisfies the initial condition is 


y= V1—-2?. 
Solution to Activity 10 


(a) The differential equation has a right-hand side 
of the form f(x)g(y), where f(x) = 1+ e” and 
gly) = 1/(1 + cosy), so it is separable. 
Separating the variables we obtain 


fo + cos y) dy = fo +e”) dz. 

Doing the integrations gives 
y+siny=x2+e"+c, 

where c is an arbitrary constant. This is the 


general solution in implicit form; it cannot be 
rearranged into the form y = y(x). 


(b) The initial condition is y(0) = 0; that is, y = 0 
when x = 0. Substituting x = 0 and y = 0 into 
the general solution, we obtain 

0+sin0=0+e° +c, 
which gives 

0=1+¢, 
so c = —1. Hence the required particular 
solution in implicit form is 


ytsiny=x+e*—1. 


Solution to Activity 11 
The given equation is y+ sin y = x + e” — 1. 
Differentiating both sides with respect to x gives 


d d 
= iny) = — Toi: 
—(y + siny) = —(e +e" — 1) 


By the chain rule, the left-hand side is 
d . dy dy 
— iny) = (1 A 
P +siny) 7” = (1+ cosy); 
The right-hand side is 


d 
— *_1)=1+¢e?". 
atte ) +e 


Hence we obtain 


dy 
1 —=1 ae 
( + cosy) + +e 


Since —7 < y < 7, we have 1+ cosy > 0, so we can 
divide through by 1 + cosy to give 

dy t1+e 

dz 1+ cosy’ 
which is the required differential equation. 


Solution to Activity 12 
The given equation is y + ln y = (1 — 2z)e?”. 
Differentiating both sides with respect to x gives 


Ly +lny) = £ ((1 — 2a)e%*) . 


By the chain rule, the left-hand side is 
d dy 1\ dy 
Sas ] — = | 1+- | —. 
g" mares ( +=) dx 

The right-hand side is, by the product rule for 
differentiation, 


a (a — 2x)e**) = (1—2z) Ge?) + (—2) (eè?) 


dz 
= (1—6z)e**. 


Hence we obtain 


(1 + 2) Ue (1 — 6x)e*. 
yj da 


Solutions to activities 


Multiplying both sides by y/(y + 1) gives 
dy — y(1— 6x)e** 


da y+1 
as required. 


? 


Solution to Activity 13 
(a) The differential equation is 


dy 3 

SS > 0). 

ie TY OY > 9) 
So the differential equation is separable and, 
after separating the variables, we obtain 


that is, 
1 — 
Dy? =% +C, 
where c is an arbitrary constant. To make y the 
subject of this equation, first take the 
reciprocals of both sides and then divide 
through by 2, to obtain 
1 
2 
Aa 2(x +c) 
Since y > 0, the general solution in explicit 


form is therefore 
1 


(e era 


The initial condition is y = 
Substituting x = 0 and y = 
solution, we obtain 

1 il 


2 V2e 
so c = 2. Hence the required particular solution 
is 


when x = 0. 


NINI = 


in the general 


1 
= aerd 


(b) We can rearrange the differential equation as 


follows: 
du B x4 
dx u?’ 


So the equation is separable and, after 
separating the variables, we obtain 


[a= fatan; 


that is, 
1,3 _1,5 
ZU =X +C, 

where c is an arbitrary constant. 
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Making u the subject of this equation, we 
obtain the explicit form of the general solution 
as 

u= (Be +0 
where b = 3c is another arbitrary constant. 


ye 


The initial condition is u = 1 when xz = 1. 
Substituting x = 1 and u = 1 in the general 
solution, we obtain 
1/3 

1= (+b), 
so b= 2. Hence the required particular solution 
is 

_ (3,5 , 2\1/8 

u=(ga°+3)". 

Solution to Activity 14 


The differential equation is 
dy _ y? cos x 


— = >0). 

ie ga 
So the equation is separable and, after separating 
the variables, we obtain 


1 
_ dy = f coszda; 
y 


that is, 


iL 1 
[(G+3) dy = f cos.zde. 
Yy y 


Doing the integrations, we obtain 


Lo g 
In|y| — 7 = sie +c, 


where c is an arbitrary constant. Since y > 0, this 
gives 


1 ! 
lny -— - = sing +c, 
y 


where c is an arbitrary constant. This is the general 
solution in implicit form, and it cannot be 
rearranged to make it explicit. 


To find the particular solution we substitute the 
initial condition, y = 2 when x = 7/2, into the 
general solution. This gives 


In2—$ =sin(7/2) +c = 1 +c, 
so c = -3 +1n2. Hence the particular solution, in 
implicit form, is 


1 
lny- = = sing — $ + ln2. 
y 
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Solution to Activity 15 


(a) 


The differential equation is 
dy 

Z =—ky/? (y>0). 
a yt (y>0) 
Separating the variables, we obtain 


pre dy = fv dt; 


that is, 

2y? = —kt +e. 
This is the general solution in implicit form, 
where c is an arbitrary constant. 


The initial condition is y = yo when t = 0. 
Substituting t = 0 and y = yo into the general 
solution, we obtain 


2y 2 =. 
Hence the required particular solution is 


2y"? = —kt + 2y, 
which can be written as 
kt = 2(yq/? — y'/?). 


The tank becomes half-empty when y = yo. 
This occurs when 


kt = 2 a = (ivo) . 


This gives 
gl 2 AP, 
kt — 2Yo == 21/2 0 > 
that is, 
,_ =v DVB 
= 7 , 


as required. 


Solution to Activity 16 


(a) 


The differential equation is 


Separating the variables gives 


1 
—dP= | Kdt. 
e 


Doing the integrations, and using the fact that 
P > 0, gives 


la P= Kt +c, 


where c is an arbitrary constant. This is the 
general solution in implicit form. 


To obtain the general solution in explicit form, 
we can apply the exponential function (the 
inverse of the natural logarithm function) to 
each side. This gives 

elnP = gare, 


that is, 

P=ektte, 
By an index law, we obtain 

P = eeft., 
Since c is an arbitrary constant, e° is also an 
arbitrary constant, and we can denote it by A. 
The value of A must be positive since e° > 0 for 
each real number c. 
Hence the general solution is 

P= Ae**, 
where A is an arbitrary positive constant. 

(b) The initial condition is P = Py when t = 0. 
Putting t = 0 and P = FP into the general 
solution, we have 

P) = Ae = A. 
Hence A = Po, and the required particular 
solution is 

P = Poe**. 


Solution to Activity 17 


(a) The initial population size is Py = 5.1 x 10°, 
and the proportionate growth rate is 
K = 0.032 year~!. The population size P after 
time t years is given by 
P = Pye*" 
that is, 
P= 1000n 
(i) After 10 years the population size will be 
P=5.1x 106e0-032x10 
=p xie 
= 7.0 x 10° (to 2 s.f.). 


That is, it will be 7.0 million, to 
two significant figures. 


Solutions to activities 


(ii) After 100 years the population size will be 
P = 5.1 x 19% e9-032 x 100 
= 5.1 x 10e? 
= 1.3 x 108 (to 2 s.f.). 
That is, it will be 130 million, to 
two significant figures. 
(b) The population size will reach 50 x 10° at the 
time ¢ years given by 
5.1 x 10ĉe0-03%t = 50 x 10°: 


that is, 
e0-032t — 50 
5.1 
Taking the natural logarithm of each side, we 
obtain 
50 
0.032t = In| — ]; 
n (2) i 
that is, 


1 50 
= in 2) = ae, 
i zar (2) (Lae 


Hence it will take 71 years (to the nearest year) 
for the population size to reach 50 million. 


Solution to Activity 18 
Suppose that the initial population size, at 
time t = 0, is Py. The proportionate growth rate is 
K = —0.03. Hence the population size P at time t 
is given by 

P = ies 
Consider a ten-year time interval from time t to 
time t+ 10. The population sizes at the start and 
end of this time interval are 

Pye79-3t and Pye 9-03(¢+10) 
respectively. 


Hence the proportion of the population size that 


will remain at the end of the time interval is 
Poe 0-03(t+10) e—0-03t ¢—0.03x 10 


Poe 00t e—0-03t 


e—0:03x 10 
— e703 
= 0.7408.... 


So a proportion of 0.74, that is, 74% (to 2 s.f.), of 
the population will remain at the end of the decade. 
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Solution to Activity 19 


For the population described in Example 9, the 
proportionate growth rate is K = 0.04 year™t. 
Hence the doubling time is 


In2 
T= — = 17. ayy 
T 7.3 years (to 3 s.f.) 


For the population described in Activity 18, the 
proportionate growth rate is K = —0.03 year™t. 
Hence the halving time is 


In2 In2 
PS. oS Si years (to 3 s.f.). 
—0.03 0.03 


Solution to Activity 20 


(a) By the equation in the box above the activity, 
the decay constant k corresponding to the 
half-life T = 2.62 hours is 


ln2 In2 
a= 0.265 hour! (to 3 s.f.). 


T 2.62 
(b) The amount m of a sample of silicon-31 that 


will remain after time t is given by 


m= moe **, 


where k is the decay constant found in part (a), 
and mp is the amount of silicon-31 at the start 
time t = 0. 


Hence the amount that will remain after 
6 hours is 


ekx6 6k 


mo = Moe 


It follows that the proportion that will remain 
after 6 hours is 


—6k 
Moe — e7 6k 
mo 
— 9 6x0.2645... 
= 0.2044.... 


So, to two significant figures, 20% of the original 
amount of silicon-31 remains after 6 hours. 


Solution to Activity 21 


The decay constant k corresponding to the half-life 
T = 5570 years is 
In2 
k = — 
T 


= 1.2444... x 1074 
= 1.24 x 1074 year~* (to 3 s.f.). 
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Let us assume that the animal died at time t = 0. 
Then the amount m of carbon-14 that will remain 
after time t is given by 


m = moe *, 


where k is the decay constant found above, and mo 
is the amount of carbon-14 at the time of the 
animal’s death. 

So the proportion of carbon-14 that will remain 
after time t is 


moe7** 


=e Kt. 


Mo 
Hence the time t at which the proportion of 
carbon-14 remaining is 85%, that is, 0.85, is given 
by 
e™™t = 0.85. 
Taking the natural logarithm of both sides gives 
—kt = ln 0.85; 


that is, 
ln 0.85 
~ K 
_ In 0.85 
~ 1.2444... x 10-4 
= 1305.971.... 


So the animal died approximately 1300 years ago. 


Solution to Activity 22 
(a) The differential equation is 
dT 


dt 
Separating the variables gives 


1 
[rr =- [rae 


Doing the integrations, and using the fact that 
T — T, > 0, we obtain 

la(T — Ta) = —At +c, 
where c is an arbitrary constant. This is the 
general solution in implicit form. 


To obtain the general solution in explicit form, 
we can apply the exponential function to each 
side. This gives 

eln(T—Ta) _ e Atte. 
From this we obtain 

T — T, = efe ™; 
that is, 

T =T, tee. 


Since c is an arbitrary constant, e° is also an 
arbitrary constant, and we can denote it by A. 
The value of A must be positive since e° > 0 for 
each real number c. 
Hence the general solution is 

T = Ta Ae 
where A is an arbitrary positive constant. 
The initial condition is T = To when t = 0. 
Putting t = 0 and T = Tọ into the general 
solution gives 

To = T, + Ae® 

=T, +A. 

Hence A = To — Ta, and the required particular 
solution is 

Ta ey Te, 
If t is very large, then e~** is very small 


(remember that A is positive), and we can 
approximate it by 0. Hence we have 


T ~% Ta + (To — Ta) x 0 
= Ts 
Hence in the long term the temperature of the 


body is equal to the ambient temperature of the 
surrounding medium. 


(This agrees with our experience.) 


Solution to Activity 23 


(a) 


The solution to Newton’s law of cooling is 
T = Ta + (%—-Ta)e™, 
where 
e T is the temperature of the cooling body at 
time t 
e To is the temperature of the body at 
time t = 0 
e T, is the ambient temperature 
e A is a positive constant. 
We know that T, = 20. Let t = 0 be the time at 
which the pathologist arrived. Then To = 30. 


So the solution to Newton’s law of cooling 
becomes 


T = 20 + (30 — 20)e™%; 
that is 
T = 20 + 10e7~. 


(b) 


(a) 
(b) 
(c) 


Solutions to activities 


We now determine the value of A. 


At time t = 1 (that is, one hour after the 
pathologist arrived), the temperature of the 
body was T = 28.5. So 


28.5 = 20 + 10e7>*1; 


that is, 

8.5 = 10e~>. 
Hence 

0.85 =e. 


Taking the natural logarithm of both sides gives 
In0.85 = In(e~*); 

that is, 
ln 0.85 = —A. 

Hence we obtain 
A = — ln 0.85. 


So the equation for the temperature of the body 
as a function of time is 


T = 20 + 10e0™0-85)t 
= 20 + 10(el9-85)¢ 
= 20 + 10 x 0.85%, 
as required. 
To find the time at which the person died, we 
need to find the value of t that gives T = 37. 
This value of t is given by 
37 = 20 + 10 x 0.85". 
Rearranging, we obtain 
0.85* = 1.7. 
Taking the natural logarithm of both sides gives 
t1n0.85 = ln 1.7, 


So 
ln 1.7 
= = —3.265.... 
ln 0.85 
So the person died three and a quarter hours 
before the pathologist arrived, to the nearest 
quarter of an hour. 


Solution to Activity 24 


The equation dy/dx + «°y = z° is linear, with 
g(x) = x? and h(x) Ss". 


The equation dy/dz = zv sin x is linear, with 
g(x) = 0 and h(x) = zsin z. 

The equation dz/dt = —3z!/? is not linear, 
because of the z!/? term. 
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(d) The equation dy/dt + y? = t is not linear, 
because of the y? term. 

(e) The equation x (dy/dx) + y = y? is not linear, 
because of the y? term. 

(£) The equation (1 + x?)(dy/da) + 2ry = 32? is 
linear, since we can divide through by 1 + x? to 


obtain 
dy 2x 3r? 
qo oe | = EE 
da l+z2 1l+<2 
which is of the required form with 
2x 3r? 


Solution to Activity 25 


(a) The differential equation has the same form as 
equation (20), with A = 1 and h(x) = e?”. So 
the general solution is given by equation (22): 


= 0 ( f eas] 


= e= (te 40) 
i2 = 
=e" ee T, 
where c is an arbitrary constant. 


(b) Differentiating the general solution we get 


dy 2 28 =T 
— = fe" — cee”. 
dx 3 

Substituting this and the general solution into 
the left-hand side of the differential equation 


gives 
dy T =p x =F 
hae a (2e — ce”) + (4e*" + ce“) 
= eE 


which is equal to the right-hand side. So 
y = že?” + ce™” satisfies the differential 
equation. 


Solution to Activity 26 


(a) The given differential equation can be 
rearranged into the form of equation (23) as 
follows: 
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The integrating factor is 


nasen (f ar) zgn 


Hence the general solution, equation (25), is 


qse (freer) 


Integration by parts gives 


y= e” (2e + f e7 ar) 


=e ( te *~—e 74 c) 


=c — zr- l, 


where c is an arbitrary constant. 


(b) After division by x (where x > 0), the 
differential equation becomes 


dy 3 =i 

de a?’ 
so 

g(x) =-- and h(x)=1. 


where c is an arbitrary constant. 


Solution to Activity 27 


(a) The differential equation can be rearranged into 
the form of equation (23) as follows: 
—-—zu=0, 


dx 


sO 
g(x) =-x and h(x) =0. 
The integrating factor is 


pa) = exp ( [ (2) ar) 


= exp(-2?/2) 


=e /2, 


Hence the general solution is 


u = e? /2 (foa) 


= ce””/ 
where c is an arbitrary constant. 
Substituting the initial condition u(0) = 2 
(that is, u = 2 when x = 0) gives 


2= ce®, 
so c= 2. Hence the solution of the initial value 
problem is 

i= 2e7”/ 2 


(b) After division by t, the given equation can be 


written as 
dy oe 2 
T a 
so it is of the form in equation (23) with 
2 
g(t) = 7 and h(t)=t. 


The integrating factor is 


p(t) = exp (J : ar) 
= exp(2lnt) (since t > 0) 
= exp (In (e) 
=i. 


Thus the general solution is 


smeed Joa) 


=G eg) 

= H + cet ?, 
where c is an arbitrary constant. 
Substituting the initial condition y(1) = 1 
(that is, y = 1 when t = 1) gives 


I= i TG 
so c= 3, Hence the solution of the initial value 
problem is 


y = qt t. 
Solution to Activity 28 


(a) The differential equation is 
dT 
— + rT = ATu e". 
dt T 0e 
Comparing this with equation (23), we see that 
the differential equation is linear with 


glt) =à and Alt)=AT ye ™. 


Solutions to activities 


Hence the integrating factor is 


p(t) = exp (J adr) Ten: 


so the general solution is 


T=e* (J AL wee ™ ar) ; 
=e (ATu few ar) ; 


(A—k)t 
= —At T CN A 
(a (a 0 X — k + ’ 


Ei AT a0 mails 
Ak 
where A is an arbitrary constant. 


(b) When k = 0, the ambient temperature at time t 
is 


4 Ae, 


T(t) = Tao eee Tao, 
and so it is a constant. Substituting k = 0 into 
the general solution we obtain 
ATao g Oai 
A — 0 
But since the constant ambient temperature in 
Activity 22 was represented by Ta, we have 


Tao=Ta, 

so the solution in this case is 
T =T, ine, 

as in Activity 22. 


T= Ae = To tAe™. 


Solution to Activity 29 
The differential equation was 
dy y-1 
dc 
and its general solution was calculated to be 


y=] Fer. 
Substituting the general solution into the right-hand 
side of the differential equation gives 

y—1 1+czr-1 

—— = —— = C 


£ z 
Differentiating the general solution gives 
dy 
—_ C, 
dz 


so the two sides of the differential equation are 
identical, as required. 
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Solution to Activity 30 


(a) This differential equation requires the 
integrating factor method. 


(b) This differential equation can be solved by 
direct integration, by separation of variables or 
by the integrating factor method. 


(c) This differential equation can be solved by 
separation of variables or by the integrating 
factor method. 


(d) This differential equation requires the 
integrating factor method. 


(e) This differential equation can be rearranged in 
the form 


dy 2 
Pee (acl Vad 
dt GDE i 


so it can be solved by the integrating factor 
method. 


(f) This differential equation can be rearranged in 
the form 


du 2 
a (u — 2)(x + 2), 


so it can be solved by separation of variables. 


Solution to Activity 31 


(a) The differential equation dr/dt = —k is directly 
integrable, so the general solution is 


r= [chat 


= —kt + c, 
where c is an arbitrary constant. 

(b) Since the initial radius of the snowball is 
100cm, the initial condition is r = 100 when 
t= 0. Substituting this into the general 
solution, we obtain c = 100. Hence the required 
particular solution is 

r = 100 — kt. 
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(c) The initial volume of the snowball (in cm?) is 


$n x 100°. If we denote the radius of the 


snowball after 2 days by R, then 

$nR? = $ x Ír x 100°, 
from which we find that 

R = (4)'/9 x 100 = 100 x 2-1/8. 
Substituting the fact that r = R = 100 x 271/3 


when t = 2 into the particular solution found in 
part (b), we obtain 


100 x 271/3 = 100 — 2k. 
It follows that 
k = 4} (100 — 100 x 27") = 50 (1 — 27"). 


(d) The snowball disappears when r = 0; that is, 
when 100 — kt = 0. This occurs when 


100 2 


It therefore takes 10 days, to the nearest day, 
for the snowball to disappear. 


Solution to Activity 32 
(a) The differential equation 


— = S 
7 CO 


is directly integrable, so the solution is given by 


p= [ teostat 


Integration by parts gives 


[ teostat =tsint — f sintat 


=tsint+cost+c, 


where c is an arbitrary constant. Hence the 
general solution is 


x =tsint+cost+ce. 
(b) The initial condition is x = 2 when t = 0. 


Substituting these values into the general 
solution gives 


2 = Osin 0 + cos0 + c= 1 +c, 
so c = 1. Hence the particular solution is 
x = tsint + cost + 1. 


Solution to Activity 33 

The differential equation is 
dg 
dt 

Separating the variables gives 


1 
[cara ftat 
x 


that is, since x > 0, 

Ing = st? FO 
where c is an arbitrary constant. To make x the 
subject of this equation, recall that exp(ln x) = z. 
Hence, taking exponentials of both sides gives 
(t? /2)+c 


=i {2 > 0): 


=g 
_ eet /2 _ Ae’ /?. 
where A = e°. The constant A must be positive, 
since ef > 0 for any value of c, but is otherwise 
arbitrary. The general solution for x > 0, in explicit 
form, is therefore 
a = Aet/ = 
where A is an arbitrary positive constant. 
The initial condition is x = 2 when t = 0. Putting 
t = 0 and x = 2 into the general solution gives 
2= Ae®, 
so A = 2. Hence the required particular solution is 
a = 2e"/?, 


Solution to Activity 34 


(a) The differential equation 
dx t 
dt cosa +a 
is separable. Separating the variables gives 


[loos +a)de = ftar 


After integration, we obtain 
sin x + $x" = st? +c, 
where c is an arbitrary constant. This is the 


general solution in implicit form. It cannot be 
rearranged to give an explicit solution. 


The initial condition is x = 0 when t = 1. 
Substituting t = 1 and x = 0 into the general 
solution gives 

0= - +c, 


=i 
sot =z. 


Solutions to activities 


Hence the required particular solution, in 
implicit form, is 
142_1 


sing + $427 = 5: 


2 
(b) To verify the particular solution, we 
differentiate both sides with respect to t: 


d,. d 
qy sine + $27) = qt — 4). 
Using the chain rule, the left-hand side is 
d d d 
qg sine + $27) = = (cosx + A 
and the right-hand side is 
d 1,2 
me — a) =t 
Thus 
(cos x + x) Ts t 
dt ` 
This equation can be rearranged to give 
ot _ @>0), 


dt cosz Fr 
which is the original differential equation, as 
required. 


Solution to Activity 35 


After dividing both sides by e% and rearranging, 
the given differential equation becomes 


dy —3t 
ae + =e š 
at 
which is a linear first-order differential equation 
with 
g(t) =1 and A(t)=e*. 


The integrating factor is 


p(t) = exp (J 1 ar) = exp(t) =<". 


Hence the general solution is 


y= e ( f tetat) =e ( f ea). 


Integrating gives 
y= e> (de® +0) 


et. E 
=ce "3e ""; 


where c is an arbitrary constant. 
Substituting the initial condition y(0) = 3 (that is, 
y = 3 when t = 0) into the general solution gives 


= AZO! . DAO A 1 
3 = ce Je = C— 3; 
soc= T. Hence the solution of the initial value 


problem is 
y = 4(Te™* — e), 


219 


Unit 8 Differential equations 


Solution to Activity 36 (b) To apply the al condition y(0) = 4, we 


(a) 
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substitute y = 7 1 and x = 0 into the sees 


First rearrange the differential equation to get : — 4 
solution. This gives 


it in the form of equation (23), 


dy n 2a 2 = a 
— sy. 
= a ’ Mi: 7 Hence the particular solution is 
= =. Hen rticular solution i 
This is a linear first-order differential equation E a : - par 
with pe x“ + 22° +1 
22 ~  A(e241) ° 
gx) = z? +1 and. lz) = which can be further simplified as follows: 
The integrating factor is therefore E (a? + 1)? T E, 
(x) [a Tape a oe 
= ex — —- dr}. 
PLL P aai 


Notice that 2x/(x? + 1) is of the form 
J (eiiie i so we can use the standard integral 


— der = In| f(2)| +c 


to Aa 
p(x) = exp(ln |z? + 1|) = exp(ln(£z? + 1)) = z? +1. 
Here we used the fact that z? +1 > 0 for all z. 


So the general solution of the differential 
equation is 


D (Je +a); 


Integrating gives 


v= (faae) 


= + ; (427+ $27? +c), 
where c is an arbitrary constant. This can be 
simplified to give the general solution 
xt + 2x? + 4e 
A(z? +1) ’ 
where c is an arbitrary constant. 


y= 
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